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ABSTRACT 

We present extensive early photometric (ultraviolet through near-infrared) and spectroscopic (opti- 
cal and near-infrared) data on supernova (SN) 2008D as well as X-ray data analysis on the associated 
Swift X-ray transient (XRT) 080109. Our data span a time range of 5 hours before the detection of the 
X-ray transient to 150 days after its detection, and detailed analysis allowed us to derive constraints 
on the nature of the SN and its progenitor; throughout we draw comparisons with results presented 
in the literature and find several key aspects that differ. We show that the X-ray spectrum of XRT 
080109 can be fit equally well by an absorbed power law or a superposition of about equal parts of 
both power law and blackbody. Our data first established that SN 2008D is a spectroscopically normal 
SN lb (i.e., showing conspicuous He lines), and show that SN 2008D had a relatively long rise time of 
18 days and a modest optical peak luminosity. The early-time light curves of the SN are dominated 
by a cooling stellar envelope (for At w 0.1 — 4 day, most pronounced in the blue bands) followed by 
^^Ni decay. We construct a reliable measurement of the bolometric output for this stripped-envelope 
SN, and, combined with estimates of i?K and M^j from the literature, estimate the stellar radius i?* 
of its probable Wolf-Rayet progenitor. According to the model of Waxman et al. and of Chevalier 
& Fransson, we derive R^^'^ = 1.2 ± 0.7 R© and i?^^"^ = 12 ± 7 R©, respectively; the latter being 
more in line with typical WN stars. Spectra obtained at 3 and 4 months after maximum light show 
double-peaked oxygen lines that we associate with departures from spherical symmetry, as has been 
suggested for the inner ejecta of a number of SN lb cores. 

Subject headings: galaxies: distances and redshifts — supernovae: general — supernovae: individual 
(SN 2008D) — galaxies: individual (NGC 2770) 
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1. INTRODUCTION 

The connection between long-duration gamma-ray 
bursts (GRBs) and broad-lined Type Ic supernovae (SNe 
Ic-bl) has recent ly piqued int erest in certain SNe as GRB 
progenitors (see lWoosle y fc Bloom 2006 for a review). A 
few SNe connected with GRBs have been fairly well ob- 
served, but there exists little early-time and high-energy 
data on the comparison set of core-collapse SNe. 

Massive stars (initial mass M > 8 M©) die vio- 
lently, and the ensuing core-collapse supernovae (SNe) 
are classified into different spectroscopic SN subtypes. 
The classification depends on the presence or absence 
of hydrogen and helium in the SN spectra, and its se- 
quence is set by the amount and kind of outer envelopes 
the progenitor retained before explosion: Type II SNe 
show prominent H, while Type lib SNe show He and 
weaker H, Type lb S Ne lack obvious H and conspicu- 
ously show He (but see lBranch et al.ll2006l ). and Type Ic 
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SNe lack both H and obvious He (see IFilippenkol Il997l 
for a review). SNe lib, l b, and Ic are collectivel y called 
"stripped-envelope SNe" (jClocchiatti et al.lll996D . While 
the number of well-observed, normal, s tripped-envelope 
SNe has grow n in the past decade (see 'Mat heson et al.l 
l2001l : lRichardson et al.lf2006:,Modiaa.2007. and references 
therein), there are few observations at very early times 
after the collapse, esp ecially when the shock breaks out of 
the stellar envelope (ICampana et ani2006t iGezari et al.l 
l200g: Sch awinski et al.ll20 08D. 

Recently, the Swift satellite (jGehrels et al.ll2004D with 
its rapid slewing capabilities, space-based nature, and 
X-ray and optical/ultraviolet telescopes has opened 
a new window for SN observations soon after their 
massive progenitors explode. A marvelous opportu- 
nity to closely study the shock breakout and evolu- 
tion of a nearby SN was offered by X-ray transient 
(XRT) 080109/SN 2008D. Durin g the Swift f ollow- 
up observations of SN 2007uy (|Nakano et al.l 120081 : 
iBlondin et all |2008c) in N GC 2770, XRT 80109 was 
discovered serendip itousl v (iBerger fc SodeTbe rg 2008; 
iKong fc Maccaronell2008l : ISoderberg et al.ll2008b[) in the 
same galaxy in data obtained with the Swift X-ray Tele- 
scope (Swift/XRT), starting on 2008 Jan. 9.56 (UT dates 
are used throughout this paper). 

After the announcement of the XRT detection, nu- 
merous g roups obta i ned data on the optical cou nter- 
part (seelPage et al| [2Q08l ISoderberg et al.ll2008bl and 
IMalesani et al.l l2009l for a detailed acc ount), which re- 
vealed it to be a SN. From spectra, I Soderberg et al.l 
(|2008aD reported a featureless continuum with a sug- 
gestion of a broad bui np near 5500 A possib ly indicat- 
ing a SN feature, and IMalesani et al.l (|2008l ) originally 
classif ied the optical transie nt as a SN Ib/c via its lack 
of H; 'Blondin etiD (12008a') and IValenti et all ()2008bf ) 
then classified the transient as a SN Ic po ssibly with 
broad spectral lines. Then we (jModjaz et al.l l2008a) were 
the first to suggest a classification as a spectroscopi- 
cally ordinary SN lb, a classification s u pported by sub- 
sequent reports (jV alenti et al. 2008a; Soderber g et al.l 
I2008bt iMazzah et a l. 2008; Malesani et al. 2009). We 
promptly established an extensive monitoring program 
of SN 2008D with facilities at Mt. Hopkins and at Lick 
Observatory, extending the observing campaign that was 
already underway for SN 2007uy in the same galaxy. 
Spectra and photometry were obtained on a nightly basis 
(weather permitting) starting 2008 Jan. 10, sometimes 
with multiple observations per night. In addition, we re- 
duced and analyzed public Swift/XRT and Swift/VYOT 
data, as well as Chandra X-ray observations. Further- 
more, we obtained optical spectra at the Apache Point 
Observatory through the awarding of Director's Discre- 
tionary Time. 

The mechanism for removing progressively increas- 
ing amounts of the H and He layers in the pro- 
genitors of stripped-envelope SNe is not fully un- 
derstood. The origin could either lie in strong 

winds of the very rna ssive progenit o r (> 30 M0; 

iChiosi fc Maedal Il986l: IWooslev et all [l993h . sudden 
eruptions ( Smith fc Owockil |2006|). interactions with 
a binary companion sta r (e.g., iNomoto et al.l 119951 : 
iPodsiadlowski et all I2004D . or the interplay of all the 
above. A larger set of SN data is needed to uncover 



trends between the underlying properties of progeni- 
tor radii, ej ecta mass, and amount of ^^Ni produced 
(|Richardsonll2009D . 

In |j2|we present X-ray data obtained with Swift/XBT 
and Chandra of XRT 080109 associated with SN 2008D, 
and in S|3| and Sj5] we present our extensive optical 
and near-infrared (NIR) photometry and spectroscopy 
on SN 2008D. Moreover, we present pre-explosion op- 
tical and NIR images from follow-up observations of 
SN 2007uy that were obtained 2.8 and 4.5 hours before 
the onset of XRT 080109, yielding stringent upper limits. 
We discuss the photometric properties of SN 2008D re- 
lated to the two phases of evolution in ij4| and construct 
bolometric light curves in fJSj fJTj compares it with the 
rest of the known SNe with observed shock breakout and 
emission from the stellar envelope, and ^Hl presents nebu- 
lar spectra that reveal double-peaked oxygen lines, which 
we attribute to global core ejecta asphericities. In ^ 
we discuss in detail the X-ray properties of XRT 080109 
compared with other kinds of X-ray events. We con- 
clude with t jlOl Throughout we dr aw comparisons 
with r e sults presented in the literature (ISoderberg et al.l 
2008bt iMaz zali et all I200I: IChevalier fc FranssonI l20'08r 



Malesani et al. 200^ and find several key aspects that 



differ. Unique aspects of our work include the very early- 
time (< 1 day after outburst) NIR data on SN 2008D, 
NIR spectra, and late-time optical spectra. 

2. X-RAY PHOTOMETRY AND SPECTROSCOPY 

Here, we discuss our independent reduction and anal- 
ysis of public Swift data and o f Chandra data (PI Poo- 
ley; iPoolev fc Soderberell200l of XRT 080109. While 
separate reductions and analysis of the same raw data 
exist in the literature (Soderb erg et al. .2Q0_8b; Li 200^; 
IMazzah et all 120081: IXu et arri2008D . our analysis shows 
that the data are consistent with a range of combina- 
tions for the contributions of the blackbody and power- 
law contributions. Constraining the spectral shape of 
the X-ray emission and its components has ramifications 
for interpreting the X-ray emission mechanism and di- 
rectly impacts the question of whet her XRT 080109 was 
produced by a jet, as for one group (jMazzali et al.|[2008l ) 
the amount of blackbody contribution and thus the small 
derived value for blackbody X-ray emitting area is a cor- 
nerstone for their jet interpretation of the origin of the 
X-ray emission. 

We downloaded the Swift/BAT and Swift/XRT 
([Burrows et al.ll2005( l data from the Swift Archive^^. The 
XRT data were processed with version 0.11.4 of the 
xrtpipeline reduction script from the HEAsoft 6.3. 
software release. We employed 4 December 2007 XRT 
calibration files. Our reduction of XRT data from 
cleaned event lists output by xrtpipeline to science- 
ready light curves a nd spe ctra is described in detail by 
iButler fc Kocevskil (|2007a[ ). In particular, our custom 
IDL tools account on a frame-by-frame basis for pileup 
and source flux lost due to hot or bad pixels and also 
perform rejection of pixels contaminated by nearby field 
sources. 

In order to best account for contamination from nearby 
sources, we extract source flux in circular regions of two 

ftp://legacy.gsfc.nasa.gov/swift/data . 
^■^ http://heasarc.gsfc.nasa.gov/docs/software/lheasoft/ . 
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Fig. 1.— Observed temporal evolution of the X-ray luminosity of XRT 080109/SN 2008D (black points) in the XRT bands (0.3-10 
keV), since tQ= 2008-01-09 13:32:49, including the Chandra observation from 2008 Jan. 19. Depending on when one defines the time of 
onset, a range of power-law slopes (green line) can be fit with a = —3.4 ± 0.6 (where Lx oc t") for data after ~ 300 s. We adopt the 
count rate to luminosity conversion from the power-law fit to the spectrum, which is 4.3 X 10^^ erg per count s~^ . The epochs on 
which optical and NIR data were obtain ed are indicated. See text for details. F or comparison, we also plot the Swift X-ray light curve 
of XRF 060218/SN 2006aj (blue points. [CTmpana et al"]|2006l : IButler et al.ll200^ . which Swift/XRT started observing ~ 100 s after the 
gamma-ray trigger. 



sizes: 16 and 4.4 pixels radius, corresponding to 90% and 
50% source containment, respectively. The spectra were 
fit using ISIS^^. Quoted uncertainties below correspond 
to 90% confidence intervals. 

2.1. X-ray Photometry: A 520 s Transient 

The X-ray transient ijBerger fc Soderber^ |2008( ) rises 
in a time 50 ± 30 s, and has a duration Tgp = 470 zb 30 s 
or A^FW HM = 80±30 s. Consistent with Soderberg et al.l 
(|2008bD . we adopt ta^ 2008-01-09 13:32:49 (± 5 sec) as 
the time of onset of the XRT — i.e., the start of the 
Swift/XBT observations. 

In Figure [1] we plot the X-ray light curve of 
XRT 080109/SN 2008D, after adequately subtracting the 
other three X-ray field sources as resolv ed by Chandra ob- 
servations (jPoolev fc Soderbergj 120081 : see below). Prior 
to 1 ks but after the peak time, the X-ray luminosity (Lx) 
decay can be fit by a broken power law in time with late- 
time index a = —3.4 ± 0.6 (where Lx oc t"), while data 
taken At < 300 s are best described by a decaying power 
law with a — —0.8 ±0.2. For comparison, we plot in Fig- 
ure [1] the X-ray hght curve of XRF 060218/SN 2006ai, 
the most recent GRB/XRF-SN (Camp ana et al.l | 200f 
Mirabal et aLll2006l:lModiaz et ah 2006: Plan et aLlbooe 



Sollerman et al. 1 120061) ■ which was discov ered by Swift 
(|Campana et al.l 120061 : IButler et all l2006f). W hile the 
expression "X-ray flash" (XRF, see lHeise et al.ll200l|) is 
used in the literature when referring to the high-energy 
event 060218, the transient 060218 was qualitatively dif- 
ferent from classical X-ray flashes, in that it likely exhib- 
ited a shock breakou t component (jCampana et al.ll2006l : 
IWaxman et al.|[2007l ). In §|9]we discuss in detail the prop- 
erties of the peculiar XRF 060218 and of XRT 080109 
compared to other, classical XRFs. 
The peak X-ray luminosity of XRT 080109 is - 700 

25 http://space.mit.edu/CXC/ISIS/ . 



times lower than that of XRF 060218 associated with 
SN 2006aj. The fluence of the X-ray transient is about 
3.0 X 10-8 erg cm"' in the 2-10 keV band (3.7 x 10"^ 
erg cm~^ in the 2—18 keV band), where we have adopted 
the count-rate conversion for the pure power-law spectral 
fit (see ij2.2p of 4.0 x 10~^^ erg cm"^ s"^ , which corre- 
sponds to 4.3 X lO''^ erg s"^ for the adopted distance to 
XRT 080109 of D=31 Mpc 

Chandra observation s (jPoolev fc Soderberd 120081 : 
iSoderberg et al.l l200 8br) have shown that three bright 
field sources lie near the X-ray transient and likely con- 
tribute substantial counts at late times even for small 
Swift/XBT extraction regions. These sources (XI, X2, 
and X3; see lPoolev fc Soderberg||200l are 3'.'8, 1'.'4, and 
2'.'6 (9.0, 3.2, and 6.1 XRT pixels), respectively, from 
the X-ray transient. Using the Chandra spectral fit- 
ting results and the Swift/XBT point-spread function 
(PSF) from the calibration database, we determine a field 
source contamination count rate of 2 x 10^'^ (6.9 x 10^"*) 
counts s~^ for the 16 (4.4) pixel extraction region. 

Correcting to the true source flux level using the PSF 
model, this shows that the flux from XRT 080109 (Fig.[l]) 
is dominated by the field source flux for any observations 
that start 5 ksec after to (i.e., the data point at 16 ksec) 
and is likely completely overwhelmed by the field source 
flux after 1 day. We therefore use the 16 pixel extraction 
region for the temporal/spectral analysis prior to 5 ksec, 
while we use the 4.4 pixel extraction region to measure 
or limit the flux at later times. 

We determine the maximum-likelihood X-ray centroid 
position and frame offset relative to SDSS (Aa =0.8" 
and AS=2A ", using 5 X-ray /optical source matches) as 
described bv lButleil (|2007bf) . Our best-fit XRT position 
is a = 09''09'"30.78% S = -h33°08'20.1" (± 1.6", J2000, 
90% conf.), which is more uncertain than the optical 
position of SN 2008D with position end figures 30.65^ 
and 20.3". These positions are consistent with those in 



4 



ISoderberg et al.l (|2008bf ). 

2.2. X-ray Spectral Fitting: Power Law vs. Blackbody 

We restrict our spectral fitting to the time period 
0-520 s after the onset of XRT 080109, where the 
contamination from X-ra. y field sources is negligible 
(|Poolev fc Soderberg|[200l . 

Different groups, performing analysis of the same 
Swift /XBT dat a, com e to different conclusions: 
ISoderberg et al.l (|2008b[ ) exclude an absorbed sin- 
gle blackbody (BB) fit on statistical grounds (with 
X^/i^ = 26.0/17) and thus favo r a si mple power-law 
(PL) fit (with x'^/v = 7.5/17): lO ( |2008D argues that the 
data can be fit by an absorbed two-component BB fit 
(with x^/'^ — 10.4/16) as well as by a pure PL fit, but 
excludes the two-component BB fit based on physical 
grounds arguing that the inferred BB radius (5 x 10^ 
cm) is much s n ialler than the possible progenitor radius. 
iMazzali et al.l (|2008D claim that for any combined BB 
and PL fit, the BB contribution has to be small, namely 
at most 14% (whose fit gives x^/'^ = 20.6/21), and that 
the implied small BB radius {^-^ 10^° cm) is evidence for 
a GRB jet (see also SJH)- However, we find that the low 
S/N data of XRT 080109 alone do not well constrain the 
contributions from a soft, therma l excess and a harder , 
PL continuum. In contrast to IMazzali et al.l (|2008D . 
we find excellent fits {x'^/v = 12.57/23, see Figure [H) 
for a combined BB and PL fit assuming a factor unity 
ratio of PL to bolometric BB flux (i.e., 50% bolometric 
contribution from a BB, see below) without increasing 
Nh relative to that in the pure PL fit (x'^/u = 12.73/24, 
Nh = 5.2tl\l X 10^1 cm-^'). This H-equivalent column 
density Nh, which is greater than that expected from 
the Galaxy (1.7 x lO^" cm""; Dickey fc Lockman, 19901 
is consistent with the lar ge reddening value inf erred 
from the optical-NIR data (jSoderberg et al.ll2008E ® , 
given standard co n versio ns between Nh and reddening 
(jPredehl fc Schmitt 1995). 

For the combined BB-PL fit, the time- integrated PL 
photon index is P ^.itol (where N{E) oc -B"^), the 
BB temperature is kT = 0.10 ± 0.01 keV and R^^ = 
10^^ cm, i.e. a factor of 10 larger than computed by 
IMazzali et al.l (2008). For comparison, the pure BB fit 
gives xV^ = 28.56/25 with i?^^ = (1.0 ± 0.2) x 10^ cm 
and kT = 0.75±0.07keV (Table|2|). Hence, the apparent 
low i?g3 from the pure BB fit alone may be an artifact of 
incorrectly fitting a low signal-to-noise ratio (S/N) com- 
posite BB plus PL spectrum with just a pure BB model. 
For completeness, we list in Table [2] our pure BB and 
PL fits, whose parameters are broadly consistent with 
the analysis performed by the groups above; specifically, 
our pure PL fit gives a time-integrated photon index of 

p _ 9 -1+0.3 

i — Z.i_Q4. 

Although the exact fraction of soft BB emission in 
XRT 080109 is poorly constrained by the data, a case can 
be made for the possible presence of BB emission through 
comparison to XRF 060218. The high-S/N X-ray spec- 
trum of XRF 060218 during the proposed shock break- 
out phase at t < 6 ks was predominantly nonthermal 
but also requi red at high confidence a soft (kT = 0.1- 
0.2 keV) BB (ICampana et al.l 12006': Bu tler et all 120061 ) 
with a ratio of ab out unity between PL and BB flux 
(|Butler et al.ll200"6l ). If we downsample the spectrum of 



XRF 060218 for i < 6 ks to contain 380 counts as in 
XRT 080109, we find that the XRF 060218 spectrum 
is fit in a strikingly similar fashion to the XRT 080109 
spectrum: a BB fit has kT = 0.72^0^ keV and local 
Nh < 3.6 x 10^1 cm-2 {x'^/iy = 35.22/23), while a PL fit 
has r = 1.7 ± 0.3 and local Nh = 2.2li J x lO^^ cm'^ 
jv = 9.25/23). From the low x^ /v, we see that the 
PL model with Nh in excess of Galactic also effectively 
overfits the soft spectral excess in XRF 060218, which 
requires a PL plus BB given all the counts. We also note 
that the light curve of XRT 080109 is quite similar to 
that observed for XRF 060218/SN 2006aj, if we allow a 
time stretch by a factor of ^ 0.1 (Fig.[T]). We discuss the 
implications and interpretation of the X-ray data in § [9l 

3. OPTICAL PHOTOMETRY 

Optical photometry was obtained w ith the Katzman 
Auto matic Imaging Telescope (KAIT; iFilippenko et al.l 
l2001h and the 1-m Nickel telescope, both at Lick 
Observatory, and the 1.2-m telescope at the Fred 
Lawrence Whipple Observatory (FLWO). We note that 
the host galaxy of SN 2008D, NGC 2770, was moni- 
tored by KAIT (|Li fc Filippenkol l2008l ) during the reg- 
ular co urse of the optical Lic k Observatory Supernova 
Search (jFilippenko et al.ll2001l : [Fihppenko 2005) as part 
of the nearby target galaxy sample, and that the regu- 
lar SN detection software successfully detected this SN 
in data from 2008 J anuary 11.42, which is 1.85 day after 
the X-ray outburst (|Li fc Filippenkoll2008D . 

3.1. Photometric Calibrations 

For photometric calibrations, the field of SN 2008D 
was observed in UBVRI during photometric nights on 
2008 Jan. 18, 19, and 20 with KAIT, in BVRI on 2008 
Jan. 12 with the Nickel 1-m telescope, and in BVr'i' on 
2008 Ja n. 18 with the FLWO 1.2-m telescope. About 
a dozen Landolt' (199^ standard-star fields were ob- 
served at different airmasses throughout each photomet- 
ric night. Photometric solutions to the Landolt standard 
stars yield a scatter of ~ 0.02 mag for all the filters for the 
Nickel telescope, and about 0.03 mag for KAIT. The SN 
2008D field was also observed for several sets of UBVRI 
images with different depth in the photometric nights. 
The photometric solutions are used to calibrate a set of 
local standard stars in the SN 2008D field as listed in 
Table |3l and a finder chart is shown in Figure [3l The 
calibration based on the FLWO 1.2-m telescope obser- 
vations ("FLWO calibration") is consistent within the 
uncertainties with the Lick calibration (from the KAIT 
and Nickel telescopes) in the filters that they have in 
common {BV). 

3.2. Lick Observatory Data Reduction 

KAIT followed SN 2008D in BVRI and unfiltered 
mode nightly (weather permitting) after its discovery. 
As SN 2008D occurred on a spiral arm of NGC 2770, 
we use image subtraction to remove the contamination 
of the host-galaxy emission. For template subtraction, 
we used pre-explosion images obtained with the FLWO 
1.2-m telescope during the course of follow-up photome- 
try for SN 2007uy, which also occurred in NGC 2770 (as 
shown in Figure [3]). 
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Fig. 2. — Observed X-ray spectrum of XRT 080109 (crosses) in the Swift/XKT bands (0.3—10 keV) averaged over the 520 s duration of the 
event. Top: Best-fit blackbody (BB; black solid line), power-law (PL; dashed black line) and combined BB-PL (with equal contributions 
- red solid line) models. The pure PL and the combined PL-BB models are both excellent fits and are nearly indistinguishable. The 
combined PL-BB fit has different parameters than the individual pure PL and BB fits. Bottom: Data residuals (crosses) from the BB 
model fit (solid line), from the PL model fit (dashed line), and from the combined BB-PL fit (red solid line). The value for the deviation 
of the respective model fit from the data is the distance between the crosses and the respective lines. 



For KAIT, the image subtraction and subsequent pho- 
tometry reduction are performed with the KAIT pho- 
tometry pipeHne (M. Ganeshahngam et al., in prepa- 
ration). Two independent image-subtraction routines 
were employed in the pipeUne. One routine is based 
on the ISIS package (jAlard fc LuptonI [l998i) as modi- 
fied by Brian Schmidt f or the High-z Supernova Search 
Team ([Riess et al.lll998[ ). a nd the other i s based on the 
IRAF26 task PSFMATCH (iPhillips &: Da vis 1995). PSF 
fitting photometry is performed on the subtracted im- 
ages, and the results from the two routines are averaged 
whenever appropriate. Artificial stars are injected into 
the original KAIT images and extracted from the sub- 
tracted images to estimate the scatter of the measure- 
ments. To convert the KAIT instrumental magnitudes 
into the standard Johnson BV and Cousins RI system, 
the color terms for the KAIT filters as determined from 
many photometric cali brations are u sed i n the conversio n 
(as detailed by Modia z et aLll2001l and iLi et al] l2001f ). 
The final error bars for the magnitudes are the scatter 
in the artificial-star simulation and the calibration error 
added in quadrature. 

We note that the FLWO 1.2-m telescope template im- 
ages were taken with the Johnson BV and Sloan r'i' 
filters, while the KAIT observations of SN 2008D were 
taken with the Johnson BV and Cousins RI filters. The 
difference between the r'i' and RI filter transmission 
curves raises the possibility of a systematic uncertainty 
during the image-subtraction process. We investigated 
this issue with the artificial-star simulations and con- 
cluded that the systematic uncertainty is smaller than 
~ 0.03 mag (which is not reported in the final uncertain- 
ties). 

SN 2008D was remotely observed with the Lick Obser- 

IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 3. — A finder chart for the local standard stars in the field of 
SN 2008D showing both SN 2008D and SN 2007uy in NGC 2770. 
The field of view is 10' X 11'. North is up and east is to the left. Dis- 
played is the r'-band image taken with the FLWO 1.2-m telescope 
on 2008 Jan. 29. 



vatory 1-m Nickel telescope in BVRI on 2008 Jan. 11 
and 12. These data were processed with proper bias and 
flat-field images, and the image subtraction and photom- 
etry were performed with the KAIT photometry pipeline 
(but modified to deal with the Nickel telescope images). 
The FLWO 1.2-m telescope template images were used to 
generate the subtracted images, and the Lick calibration 
and the proper color terms for the Nickel telescope fil- 
ters are used to convert the instrumental magnitudes into 
the standard system. The final photometry of SN 2008D 
from the KAIT and Nickel observations is listed in Ta- 
ble |2 
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3.3. FLWO 1.2-m Telescope Data Reduction 

The field of SN 2007uy in NGC 2770 was monitored 
by the FLWO 1.2-m telescope in the context of the CfA 
SN monitoring^^ efforts when SN 2008D was discovered, 
providing us with valuable pre-explosion images. We fol- 
lowed SN 2008D on a nightly basis (weather permitting) 
in BVr'i' . For performing photometry of FLWO 1.2-m 
telescope images, we adopted the image- analysis pipeline 
of the SuperMACHO and ESSENC E col laborations (see 
Rest et all l2005l iGarg et all l2007l and iMiknaitis et all 
20071 for details) . The operations of this optical pho- 
tomet ry pipeline are discussed in depth bv lHicken et al.l 
(|2009) . In brief, we employed differential photometry by 
measuring the brightness of the SN with respect to a set 
of stars (see Figure[3l) in the SN field and we employed the 
DoPHOT photometry package (jSchechter et al.lll993D to 
measure the flux of the SN and its comparison stars. We 
performed image subtraction with the robust algorithm 
of Alard & Lupton (1998; see also^l_ard 2000, using as 
templates the 1.2-m images taken on 2008 Jan. 08 of 
the field. The "FLWO calibr ation" and the color terms 
for Keplercam (|Modiajl2007f l were used to convert the 
instrumental magnitudes into the standard system. 

We present the 1.2-m BVr'i' photometry in Tabled 
For multiple observations on the same night we report the 
weighted mean of those observations. Intranight obser- 
vations did not reveal statistically significant variations 
for SN 2008D. We furthermore derive upper limits for 
SN 2008D based on 2008 Jan. 9.44 images of the field of 
SN 2007uy, three hours before the onset of XRT 080109 
(§[2]), and present them in Tabled The upper limits were 
estimated by placing fake stars with a range of magni- 
tudes at the position of SN 2008D on the template images 
and then determining when they are detected at varying 
thresholds above background. 

3.4. UVOT Data Reduction 

The Swift/ VYOT has followed SN 20 8D intensively 
(iKong et a l. 2008b: Immler et al. 2008a; Li e t al.l 120081 : 
llmmler et al, .2008b; Soderberg et al. 2008bl We re- 
trieved the level-2 UVOT (jRoming et al.ir2005l ) data for 
SN 2008D from the Swift data archive in all available fil- 
ters (i.e., in the UVW2, UVM2, UVWl, U, B, and V 
filters). 

We used the UVOT observations of SN 2007uy on 2008 
Jan. 6 as the template images to reduce the data of 
SN 2008D (a total exposure time of 808 s in each of 
the filters). For the [/, S, and V filters, we employed 
th e photonietric r eduction and calibration as described 
bv lLi et al.l (|2006[ ) to perform photometry of SN 2008D 
on the subtracted images, while for the Far-UV filters 
{UVW2, UVM2, UVWl), we used a 3 " aper ture and 
employed the calibration of iPoole et al.l (|2008f ) . Many 
tests have shown that t he tw o calibrations by ' Li et al.l 
([200d ) and iPoole et all ()2008f ) are consistent with each 
other. 

The UVOT photometry of SN 2008D is reported in 
Table [5l and is presented in Figure [H along with the 
overall optical photometry, with respect to the start 
of XRT 080109 (©. The UVOT ^U values reported 
here arc fully consistent with those in ISoderberg et al.l 



See |http: / /www.cfa.harvard.edu/supernova/index.html 



(|2008bH . while our U photometry is systematically fainter 
by 0.2 mag (consistent within systematic errors, see be- 
low), and our UVWl, UVW2, and UVM2 3cr upper lim- 
its are systematically fainter by 0.2-1 mag. Our analysis 
used larger binning periods than that of Soderberg et al., 
producing deeper images and enabling additional detec- 
tions in UVWl. The S/N of the UV detections is poor, 
and thus precision photometry difficult. Specifically, the 
choice of aperture size affects the final SN magnitude; for 
example, using a 2" aperture for the UVWl band gives 
magnitudes that are ~ 0.23 mag brighter than for the 3" 
aperture that we adopt (even after appropriate aperture 
corrections). Thus, we estimate the systematic uncer- 
tainty of the final UV photometry to be on the order of 
^ 0.25 mag. 

3.5. Near-Infrared Photometry 

We obtained NIR (JHKs) photometric measurements 
with the refurbished and fully automated 1.3-m Pe- 
ters Automated Infrared Telescope (PAIRITEL)^^ lo- 
cated at FLWO. PAIRITEL is the world's largest NIR 
imaging telescope dedicated to time-domain astronomy 
(Bloom et al. 2006). The telescope and camera were 
part of the Two M icron All Sky Survey (2MASS; 
ISkrutskie et ~aIl l2006D project . H-, and iC^-band (1.2, 
1.6, and 2.2 /Ltm: ICoheii et alJ 2003) images were acquired 
simultaneously with the three NICM0S3 arrays with in- 
dividual exposure times of 7.8 s. 

The PA IRITEL reduction pipeline software 
(jBloom et al., ,2006i) is used to estimate the sky 
background from a star-masked median stack of the SN 
raw images. After sky subtraction, the pipeline is used 
to cross-correlate, stack, and subsamplc the processed 
images in order to produce the final image with an 
effective scale of 1" pixel""'^ and an effective field of 
view (FOV) of lO'xlO'. The effective exposure times 
for the final "mosaic" images ranged from 15 to 20 min. 
Multi-epoch observations over the course of the night 
(totaling between 40 min to 4 hours) were obtained to 
search for intranight variations, but none were found at 
a statistically significant level. 

For the NIR photometry, we used the image-analysis 
pipehne of the ESSE NCE and SuperMACHO p rojects 
ilRest et all [20(351: IGarg et all l2007l : IMiknaitis et all 
2007t). We performed d ifference photometry (following 
Wood-Vasev et al.ll2008h using the image-subtraction al- 
gorithm of I Alard fc LuptonI (|l998f) . similar to our opti- 
cal photometry. As a template image for the field of 
SN 2008D we used our PAIRITEL images of SN 2007uy 
taken on 2008 Jan. 9. Thus, we also derive NIR up- 
per limits for SN 2008D based on images of the field of 
SN 2007uy, taken on 2008 Jan. 9 08:26:19, 4.5 hours be- 
fore to, the onset of XRT 080109 (see and present 
them in Table [6l The upper limits were estimated in 
the same fashion as for the optical data. We performed 
the final calibration onto the 2MASS system by using as 
reference stars the field stars from the 2MASS catalog, 
of which there were 20—30 in the FOV. No color-term 
corrections were required since our natural system is al- 
ready in the 2MASS system. We present our PAIRITEL 
photometry in Table [H] 

See |http://www.pairitel.org/| . 
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4. EARLY-TIME BROAD-BAND PHOTOMETRY 

Figure H shows the UV {UVW2,UVM2,UVW1,U), 
optical {BVr'i'), and NIR (JHKs) photometric evo- 
lution of SN 2008D, starting as early as 0.1 day af- 
ter burst, and going out to 3 5 day. The data revea l 
a two-component light curve (|Soderberg et al.1 l2008bD : 
the first due to the cooling stellar envelope that is most 
expressed in the blue bands (At — 0.1—4 day) and 
the second due to the standard radioactive decay of 
56Ni and ^^Co (At > 5 day) with half-hfe times of 
6.1 day and 77.3 day, respectively, which lead to en- 
ergy deposition behind the SN photosphere. A num- 
ber of theoretical studies have computed shock break- 
out in the context of SN II 1987A and correspond- 
ing optical li ght curves due to the cooling stellar en- 
fe.g., IWoosley et al.l 119871: lE nsman fc Burrowi 
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[amngiiiESEi 



1993 J (e.g., IBlin nikov et a 



Ic 2006aj (iWaxman et al 



Blinnikov et al. 2000), SN lib 
]Tl99|^nd broad-lined SN 
I2007t ). In the following 



we wi ll concentrate on the models bv IWaxman et all 
(|2007| ) as used by ISoderberg et al.l (|2008bl ) to fit to 
the Swift/\JYOT data and their ground-based data, as 
well as the most recent model bv lChevalier fc FranssonI 
(HqM)- In EH we wiU briefly discuss the implications 
of the differences between those two models for deriving 
the progenitor radius. 

The most remarkable aspects are the early-time 
Swift/WOT data and our very early optical and NIR 
data points (at At — 0.84 and 0.71 day, respectively), of 
which the NIR data are the earliest broad-band observa- 
tions reported of a SN lb. They were obtained as part 
of our observing campaign of SN lb 2007uy in the same 
galaxy. Here all four {BVr'i') filters indicate a first peak 
at 0.8-1.8 day. We note that these earliest BVr'i' data 
points are consistent with the cooling stellar envelope 
fits (dashed lines in Fig. |4]) by ISoderberg et all ()2008l3l ) 
to their ground-based data, while the fit in the UVM2 
passband does not agree with the photometry as per- 
formed by us; the fit predicts brighter magnitudes than 
observed, e.g., by ^ 0.2 mag for the first UVM2 upper 
limit. 

Our y-band light curve of SN 2008D peaks on 2008 
Jan. 27.9 ± 0.5 (i.e., at At = 18.3 ± 0.5 day) at appar- 
ent magnitude my — 17.33 ± 0.05 mag. The S-band 
rise time (between shock breakout and maximum light) 
is slightly shorter than in the V band, namely 16.8 ± 0.5 
day, peaking on 2008 Jan. 26.4 at tob = 18.5 ± 0.05 
mag. Our apparent peak magnitudes are consiste nt with 
other r eports within the u n certainties (ISo derbcr g et al.l 
I2008bt iMazzah et all l2008l : iMalesani et ai.ll2009if l. The 
rise times of SN 2008D are at the long end of the range 
of those observed f or stripped-envelope S Ne, and sim- 
ilar to SN 1999ex (jStritzinger et al.ll2002f ). After cor- 
recting for a Galactic extinction and host-galaxy red- 
dening of E{B - y)Host = 0.6 ± 0.1 mag (see § E]), 
these values correspond to a peak absolute magnitude of 
Mv = -17.0 ± 0.3 mag and Mb = -16.3 ± 0.4 mag for 
SN 2008D, using a distance modulus of ^ = 32.46 ± 0.15 
mag (which corresponds to a distance of D ~ 31±2 Mpc). 
This distance modulus is based on the measured helio- 
centric velocity of NGC 2770, corrected for the effect of 
the Virgo cluster and the Great A ttractor by us i ng th e 
local velocity field model given bv iMould et aH ([200(Th . 



These are modest peak luminosities compared to those 
of other known SNe lb and SNe Ic, in particular about 1 
mag less luminou s than the mean of the SN lb sample in 
Rich ardson et all (j2006) (< My -17.9 ± 1.0 mag, 
recomputed using Hq — 73 km s"^ Mpc"^ for those SN 
with known luminosity distances), but within the range 
of observed values. SN 2008D is fainter than SN 1999ex, 
a well-observed SN lb (by '^0.5 mag and 0.8 mag in V 
and B, respectively) and SN 1993J (by ^ 0.4 mag). 

5. SPECTROSCOPIC OBSERVATIONS 

Optical spectra were obtained with a variety of instru- 
ments: Gemini-South via queue-scheduled observations 
(GS-2007B-Q-2, PI Chen), the 6.5-m MMT, the 10-m 
Keck I and Keck II telescopes, the 6.5 m Clay Telescope 
of the Magellan Observatory located at Las Campanas 
Observatory (LCO), the 3-m Shane telescope at Lick Ob- 
servatory, the Astrophysical Research Consortium 3.5-m 
telescope at Apache Point Observatory (APO), and the 
FLWO 1.5-m t elescope. The spe ctrographs utilized were 
CM OS-South iH ook et al1[200l at Gemini, Blue Chan- 
nel (ISchmidt et al. 1989 ") at the MM T, LRIS (|Oke et all 
I1995D and H IRES ( Vogt et a l."1994') at Keck I, DEIMOS 



( Faber et a l. 2003) at Keck II, LDSS-3 (Mulchaey 
Gladders 2005) at LC O, the Kast Double Spectrograph 
UMiller fc Stongll993l ) at Lick, the DIS spectrograph at 
the A PO, and the FAST spectrograph (|Fabricant et all 
Il998l ) on the FLWO 1.5-m telescope. A detailed journal 
of observations is shown in Table [T] 

Almost all optical spectra were reduced and calibrated 
empl oying standard techniques (see, e.g., iModiaz et al.l 
HqoU) in IRAF and our own IDL routines for flux cali- 
bration. In addition, the HIRES spectrum was reduced 
using HIRcdux^^, an IDL package developed by one of 
us (J.X.P.) for the HIRES mosaic, and the GMOS-S and 
DEIMOS spectra were reduced using slightly different 
techniques (see lFolev et~alll2006l[200^ . 

Extractions of the low-resolution sp ectra were do ne us- 
ing the optimal weighting algorithm ljHornel[l986l ). and 
wavelength calibration was accomplished with HeNeAr 
lamps taken at the position of the targets. Small-scale 
adjustments derived from night-sky lines in the observed 
frames were also applie d. The spectra w ere either taken 
at the parallactic angle ()Filippenkolll982( ). at low airmass, 
or with an atmospheric dispersion corrector, in order to 
minimize differential light loss produced by atmospheric 
dispersion, except for the GMOS-South spectrum. Tel- 
luric lines were removed f o llowin g a p rocedure similar 
to that of IWade fc Hor^ (| 1988( 1 and iMatheson eFall 
(j2000b). The final fiux calibrations were derived from ob- 
servations of spectrophotometric standard stars of differ- 
ent colors to correct for second-order light contamination 
in the MMT and FAST spectrographs (see lModjaz et all 
I2008bl for more details) . 

Four NIR spectra were obtained on 2008 Jan. 14.45, 
21.51, 28.51 and 2008 Mar. 14.48 using the 3.0-m tele- 
scope at the NASA Infrared Telescope Facility (IRTF) 
with the Spe X medium-resolutio n spectrograph in prism 
(LRS) mode (|Ravner et al.ll2003D . The SpeX instrument 
provides single-exposure coverage of the wavelength re- 
gion 0.8-2.5 /im in either cross-dispersed or prism mode. 
The data were reduced and calibrated using a package 

|http: / /www.ucolick.org/ ~xavier/HIRedux/index.html| . 



8 



10 



O 



• 1— I 

bJD 



20 



I 30 



40 



jiiiiii 1 1 — p 

T 

T 

t 

5 

5 
5 



I lui 



"T 1 I I I I I 11 



"T 1 I I I I I I 



I — I I I I I llj r 

-•K,-10 

wur H-8 

J-6 




.... S08 
.... SOS 
.... 808 
.... SOS 
.... SOS 
.... SOS 
.... SOS 
.... SOS 
I 



:V 
:B 
:U 

:UVW1 
:UVM2 
:UVW2 

J I 



uvwi+e: 

UVM2+8 
UVW2+11 



-1.0 0.1 1.0 10.0 100.0 

Days after Outburst 

Fig. 4.— Observed optical and NIR light curves of SN 2008D after the onset (right panel) of XRT 080109 at to = 2008-01-09 13:32:49 
(± 5 s), which we adopt as the time of shock breakout. The filled circles show UVW\, UVM2, UVWl, U, B, V data from Swift/VVOT, 
the empty circles are BV data from KAIT (its RI data are not shown for sake of clarity), while the empty squares are BVr'i' data from 
the FLWO 1.2-m telescope. JHKs data (filled stars, triangles, and squares) are from PAIRITEL. Note the very early optical data points 
(At =0.84 day after shock breakout) from the FLWO 1.2-m telescope, as well as NIR data (at At =0.71 day) from PAIRITEL. Swift/VVOT 
upper limits are indicated by the arrows. We also plot the pre-explosion upper limits derived from the 1.2-m CfA (at At = —2.8 hr) and the 
PAIRITEL (at At = —4.5 hr) data (left panel). The data have not been corrected for extinction. We note that our earliest ground-based 
BVr'i' dat a points are consistent with the light-curve fits by Soderberg et al. (2008, 308; dotted lines), who use the envelope BB emission 
model from lWaxman et al] I I2007I) . 

of IDL routines specifically de signed for the redii ction of 
SpeX data (Spextool v3.2; C ushing et aI]|2004D . Cor- 
rections for telluric absorption were performed using the 
extracted spectrum of an AOV star and a speci a lly de - 
signed IDL package developed by IVacca et al.l (|2003D . 
Additional details of the data acqu i sition and reduction 
methods are given bv lMarion et al.l ()2008f) . 

5.1. High- Resolution Spectroscopy: Probing the 
Interstellar Medium 

We obtained one high-resolution spectrum with Keck 
I on 2008 Jan. 12.5 (i.e., at At = 2.95 day), over the 
observed wavelength range 3440-6300 A. We detect un- 
derlying H II region emission lines (Balmer lines and 
[O II] A3727), as well as Na I an d Ca II HfcK absorption 
lines in the host galaxy (see also ISoderberg et "al1l2008bl : 



iMalesani et"alll2009[ ). The Na I and Ca II H&K absorp- 
tion lines have identical redshifts as the centroid of the 
[O II] and Balmer emission hues, z = 0.00700 ± 0.00005. 
Since both emission and absorption lines show the same 
velocity, we consider the Na I and Ca II H&K absorption 
lines to have interstellar origin (as opposed to circum- 
stellar) and attribute the offset of ~160 km s^^ from 
the measured redshift of the nucleus of the host g alaxy 
NGC 2770 (^nuc = 0.006451±0.000087: [lMcoiraIll999) 
to the galaxy's rotation curve. Given that NGC 2770 
is an inclined, edge-on, typical spiral (Sb/Sc) galaxy 
()Soderberg et a l. 2008b; Thonc et al. 2009), rotation on 
the order of 160 km s^^ at the distance of the SN from 
the nucleus (~ 9 kpc) is plausible; the full rotation curve 
of N GC 2770 extends to ~400 km s"^ (jHaynes et al.l 
[1993). 
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Fig. 5. — Spectral evolution of SN 2008D, dereddened by E(B — V)Host = 0.6 mag and labeled with respect to date of shock breakout, 
to = 2008-01-09 13:32:49 (indicated by At), and to date of V-band maximum, 2008 Jan. 27.9 (indicated by tvmax)- Note the fleeting 
double-absorption feature around 4000 A in our early spectrum at At = 1.84 day, which is discussed in 15.31 Some spectra have been 
binned to more clearly visualize their effective S/N. Unfortunately, we have a gap in our spectroscopic coverage around maximum light. 
The characteristic optical He I lines (due to blueshifted He I AA4471, 5876, 6678, 7061) become visible starting At « 12 day or t^max ~ ~6 
day. Note that not all spectra are presented in this figure. For late-time spectra, see § |8] 



Since the absorption Hnes are saturated, we can only 
place lower limits to the column densities of Na I and 
Ca II H&K: iVNa i > 1.0 x lO^^ cm-" and A^ca ii > 
1.4 X IQ-'^'^ cm"". We measure an equivalent width (EW) 
for Na I D2 of 0.67 ± 0.04 A (at rest), and for Ca II K 
of 0.55 ± 0.04 A (at rest). 

5.2. Low- Resolution Spectroscopy 

In Figure [5l we present the sequence of early-time op- 
tical spectra after subtracting the SN recession velocity, 
dereddened by -B(B-F)Host = 0.6 ± 0.1 mag (§111). Fur- 
thermore, the spectra are labeled according to age (At) 
with respect to date of shock breakout {to = 2008-01-09 
13:32:49) and epoch (tymax), defined as days relative to 
V-hand maximum (2008 Jan. 27.9 = JD 2454493.4, see 
§ |4]) . Knowing both reference dates is crucial for inter- 
preting the spectra and their temporal evolution. 

These are amongst the earliest spectra presented of 
a normal SN lb (starting 1.70 day after shock break- 
out, which means 17.0 day before y-band maximum) 
and provide insights into the different layers of the SN 
ejecta. As the photosphere recedes with time, it reveals 
lower- velocity material since the SN ejecta are in homol- 



ogous expansion (where v (x R). Moreover, we obtained 
late-time spectra at tymax = +91 and +132 day that 
show the advent of the nebular epoch, in which the spec- 
trum is marked by strong forbidden emission lines of 
intermediate-mass elements [O I] and [Ca II] (not shown 
in Fig. [5]). We discuss the late-time spectra in §|S1 

At early times our densely time-sampled series shows 
the dereddened spectra to consist of a blue, nearly fea- 
tureless continuum (except around 4000 A), with su- 
perimposed characteristic broad SN lines not attributed 
to hydrogen dMalcsani ct al. 2008: Blondi n et ani2008at 
iValenti et all [2008bl : ISoderb erg ct al. 20081?). In § [O] 
below, we discuss the fleeting, double- absorption feature 
around 4000 A which is only seen in our spectrum at 
At= 1.84 day. 

Figure [5] shows that the optical helium lines 
(blueshifted He I AA4471, 5876, 6678, 7061) be- 
come apparent starting at At « 10-12 day (8-6 day 
before V^-band maximum) and grow stron ger over time, 
allowin g this SN to b e classified as a SN lb (iModiaz et al 
2008a: IValenti et al . 2008a- se e also ISoderberg et al . 
I2008d iMazzah et all I20d8l : iMalesani et al.l l2009[ ). 
This is in stark contrast to all SNe associated with 
GRBs that have been exclusively broad-lined SN 
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Ic (e.g., iGalama et al l 11998!; IS tanek et all 120031: 



Hiorth et alJ 
20061: iMirabal et al.l 



2omr 



Plan et al.ll2006D . 



Malcsani ct al. 2004 : iModiaz et all 
12006: iCampana et all l2006t 



5.3. The Double-Absorption Feature around 4000 A 

Perhaps the most intriguing aspect of the very early- 
time spectra of SN 2008D is the transient, prominent, 
double- absorption feature around 4000 A (Fig. [H]). It is 
only present in our second spectrum, the MMT-|-Blue 
Channel spectrum tak en on Jan 11.41 (At — 1.84 day), 
and was first no ted bvlBlondin et all (l2008aD . 

Furthermore, iMalesani et al.l ( 2009D obtained spectra 
of SN 2008D at nearly the same epoch, Jan. 11.32 
(i.e., within two hours of our Gemini and MMT spec- 
tra, see Table |7|) with wavelength range extending down 
to 3700 A, and note that their spectra show the same 
double-absorption feature. The double-absorption fea- 
ture is apparently absent from our very first spec- 
trum, the Gemini-f CMOS spectrum taken on Jan 11.26, 
but the spectrum is extremely noisy below ^ 4500 A 
due to the par ti cular s etup used. The spectrum of 
ISoderberg et al.l (|2008bD on the same and subsequent 
nights apparently does not extend sufficiently far to the 
blue to include this feature at '^4000 A (their Fig. 3). 
The double-absorption line is no longer visible in spec- 
tra taken > 1 day later (i.e.. At > 2.99 day) and the two 
components appear to blend into a single, broad (FWHM 
w 23, 000 km s~ ) absorption feature (maximum absorp- 
tion at rest frame ~ 4100 A), as clearly seen in the Jan. 
13.30 MMT-hBlue Channel spectr um (At = 3.74 d ay) in 
Figure (6] Also, spectra taken by IMalesani et al.l (j2009l ) 
do not show the feature in spectra taken after Jan. 13.07, 
i.e.. At > 3.51 day. 

The double- absorption feature appears to consist of 
two overlapping P-Cygni profiles, with rest-frame full 
width at half-maximum (FWHM) of - 11,500 km s"^ 
and ~ 11, 000 km s~^ for the bluer (maximum absorption 
at rest- frame wavelength ~ 3980 A) and redder (maxi- 
mum absorption at rest-frame wavelength ~ 4240 A) ab- 
sorptions, respectively. Both absorption components are 
thus significantly narrower than other absorption fea- 
tures in the same spectrum, whose FWHM ranges be- 
tween - 15,000km s"^ and - 30,000km s"\ 

In the following section we attempt to address the line 
identification for this feature. Amongst the published 
supernova spectra, the closest match to the "W" feature 
that we can find is in the earliest spectrum of SN 2005ap 
(|Quimbv et al.ll2007f ) from 2005 March 7 (~ 1 week be- 
fore maximum light), which shows a similar feature at 
slightly redder wavelengths (by ~200 A). Quimby et 
al. used the para meterized sup erno va spectrum synthe- 
sis code SYNOW (iFisher et al.l[l99l iBranch et al.ll200l 
to fit that spectrum with four ions. They ascribed the 
"W" feature to a combination of C III, N III, and O HI 
with a photospheric velocity of 21,000 km s^^. We down- 
loaded the spectrum of SN 2005ap from the SUSPECT 
database^°, generated SYNOW models using only those 
three ions, and were able to satisfactorily reproduce the 
"W" feature in SN 2005ap using similar fitting parame- 
ters to those of Quimby et al. We then increased what 

http://bruford.nhn.ou.edu/ suspect /indexl. html . 
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Fig. 6.— Early-time spectra of SN 2008D in the range 3500- 
4700 A, from At = 1.70 to At = 4.79 day past the X-ray outburst, 
not corrected for reddening. The transient double-absorption fea- 
ture is marked with arrows in the spectrum at At = 1.84 day. 

SYNOW calls the "photospheric velocity" (uphot) of the 
fit to match the observed wavelengths of the "W" feature 
in SN 2008D and found that ?;phot « 30,000 km s"^ pro- 
duced a good match to the "W" as well as to some other 
weaker features, as shown in Figure |7| (using a continuum 
temperature of ~13,000 K, consistent with the BB tem- 
perature as derived from our photometry data. Table |8|). 
SYNOW did not produce any other strong features that 
are associated with these ions, matching the full extent 
of the observed MMT spectrum. 

While our SYNOW fit produces a good match to 
the data, we note that SYNOW fits do not produce 
unique line identifications and the suggested identifi- 
cation here n eeds to be confirmed by detailed spec- 
tral modeli ng (Mazzali fc Lucv''l998': iBaron et al.lll999l : 
iDessart fc H ilher 2008). Furthermore, while the "W" 
feature is no longer visi ble in spectra of SN 2 008D taken 
~ 1 day later (see also IMalesani et 311120091 ) . it seemed 
to linger for ~ 1 week in the spectrum of SN 2005ap (at 
slightly lower velocities). 

At face value, an identification of the "W" feature with 
C HI, N HI, and O III at high velocities seems plausi- 
ble. We expect the early supernova spectrum to have 
high velocities based on an extrapolation of the observed 
velocities at late times back to the earliest observations. 
In addition, what SYNOW designates as the "excitation 
temperature" is high (Tgx « 17,000 K) at these early 
times^^ and many of the elements that will contribute 

The apparent incongruency between the excitation and 
continuum temperatures probably arises from the LTE (local 
thermodynamic equ ilibrium) assumption in the SYNOW code 
IIBranch et al.ll200a) . or because we underestimated the reddening, 
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to the later spectra are more highly ionized. The tran- 
sient nature of the feature might then be explainable 
because the elements involved will recombine to lower 
ionization stages as the ejecta cool (at ~15,000 K where 
most of doubly ionized elements recombine to singly ion- 
ized species). 

Alternatively, the reason for the short-lived nature of 
this "W" feature could be asphericity in the SN ejecta. 
We have evidence from late-time spectra for asphericity 
in at least the inner SN core (see § [5]) , which may be con- 
nected to some inhomogeneity in the SN ejecta density 
or chemical abundance profile at the outermost radii that 
the early-time spectra are probing. If indeed asphericity 
is involved, then existing spectral synthesis codes (such 
as SYNOW which assumes spherical symmetry) are chal- 
lenged to accommodate that aspect. ■^^ 
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Fig. 7. — The reddening-corrected early-time spectrum of 
SN 2008D (black) at At = 1.84 day (2008 Jan 11.41) and for com- 
parison, the spectrum of SN 2005ap (Quimby et al. 2007). The fit 
via the parameterized spectral synthesis code SYNOW is shown 
in red: the "W," along with some other weaker features, can be 
reproduced with a combination of O III and CIII/NIII. 

6. BOLOMETRIC LIGHT CURVES 

Many lines of evidence indicate significant host-galaxy 
extinction toward SN 2008D: the high- resolution spec- 
trum (§ 15. ip shows strong interstellar Na I D and 
Ca II lines at the recession velocity of the host galaxy 
(with EW[Na I D] = 0.7 A), while low-resolutio n spec- 
tra s how in addition diffuse i nterste llar bands (jHerbid 
|1975[) . FromlSoderberg^et al.' (' 2008bf ). we adopt E{B - 
y)Host=0 .6 mag for S N 2008D, and a Milky Way extinc- 
tion law (!Fitzpatrick'T999) with Ry = 3.1. This value 
is used to dcrcddcn all of o ur spectra and is c onsistent 
with the value adopted bv lMazzah et all (I2008D. but 0. 2 
mag smaller than that used by iMalesani "eTall (|2009l) . 
The choice of extinction law does not strongly affect our 
results. 

This value for the reddening is consistent with fitting 
an absorbed BB spectrum to the spectral energy distribu- 
tion (SED) of SN 2008D constructed from our simultane- 
ous broad-band {UVW2 — > Kg) photometry (see 16. ip . 

but this is less likely. 

The MMT spectra (including the one display- 
ing the double-absorption feature) will be made pub- 
licly available through the CfA Supernova Archi ve at 
|http : //www . cf a . harvard . edu/supernova/SMar chive ■ html | . 



In addition, comparison of spectra and B ~ V col- 
ors with those of other SNe lb at comparable epochs, 
such as SN 1999ex (jStritzinger et al.ll2002l: iHamuv et all 
I2002f ). yields E{B -V) = 0.6-0.7 mag. In particular, we 
compared SN 2008D with SN 2005hg, a normal SN lb 
that suffered very little host-galaxy extinction (Modjai 
l2007f ). Comparison of the spectra (Fig. [T^ and the 
broad-band SEDs constructed from our UBVr'i' JKHg 
photometry for both SNe indicates a host-galaxy extinc- 
tion for SN 2008D of E{B -V) ^ 0.65 mag, if both SNe 
have similar color evolution. 

6.1. Constructing Bolometric Light Curves 

The broad-band dataset (spanning 2200-24,100 A) 
presented here lends itself uniquely to computing the 
bolometric output of SN 2008D. We determine it in two 
different ways: (a) by fitting a BB to the broad-band 
photometry and deriving the equivalent BB luminos- 
ity {Lbb) by analytically computing the corresponding 
Planck function, and (b) by performing a direct integra- 
tion of the UVWlBVRr'Ii'JHKs broad-band magni- 
tudes {Luvwi^K^), including the far-UV filters UVW2 
and UVM2 in case of detections. In both cases, the 
Swift UV and our NIR data are crucial, as they provide 
a factor of at least two in wavelength coverage and lever- 
age for the SEDs. While theoretical considerations might 
justify using a simple BB approxim ation for the SN emis- 
sion shortly aft er shock breakout (jWaxman et al.l [20071 : 
IChevalier fcfe msson 200l), we do note that without de- 
tailed models of SN Ib/c atmospheres, it proves difficult 
to say how well BB fits approximate the true SN lumi- 
nosity over the full light curve of SN 2008D. We address 
this question more in the remainder of the section, when 
comparing the bolometric luminosities computed via BB 
fits and those via direct integration. 

We start by interpolating the light curves by 
Swift/\JYOT (in UVWl^UBV), FLWO 1.2 m (in r'i'), 
and PAIRITEL (in JHKg) onto a grid spaced regularly 
in time, as the observing times were different for different 
telescopes. We then sampled those light curves at inter- 
vals closest to the actual observing time of Swift/VYOT 
and the FLWO telescopes. After removing extinction 
{E{B—V)iiost = 0.6 ± 0.1 mag and assuming a the Milky 
Way extinction l aw parameterization with Ry — 3.1; 
iFitzpatrickl Il999f ). we converted our photometry into 
monochromatic fluxes at the effective wavelengths of the 
broad-band filters, using zeropoints from lFukugita et all 
lfT995l ) (for UBVr'i', since the UVOT UBV photometric 
system is clos e to the Johnson UBV with Vega defining 
the zeropoint: |Poole et al.ll2008l ) and lCohen et all (|2003h 
(for JHKg). We incl uded 0.0 5 mag of systematic error 
in the UVOT filter s (iLi et a l. 2006) and 4% conversion 
error when using .Fukugita et al.. (1995i) . 

For our first method, we perform a least-squares fit of a 
BB spectrum to the SED of each epoch and constrain the 
derived BB temperature (Tbb), BB radius {Rbb), and 
BB luminosity {Lbb) as a function of age (see Table [8] 
and Fig. [8|). The error bars include the uncertainties 
from the formal fits as well as uncertainties introduced 
by the uncertain reddening (±0.1 mag), which affect the 
derived BB fit parameters differently, depending on the 
underlying SED. 

We do not include the systematic uncertainty in dis- 
tance (of order -20%). At later times (Ai > 8-10 day) 
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Fig. 8.— Left: Temperature (top) and radius (bottom) derived from a BB fit to the optical and NIR SED of XRT 080109/SN 2008D 
(filled circles) as a function of time, assuming _D = 31 ± 2 Mpc and corrected for reddening (E(B — V)Host = 0-6 ± 0.1 mag). The 
error bars include the uncertainties as reported from the formal fits, as well as the effect of uncertain extinction. SN 1999ex (open circles, 
[Stritzinger ct al. 2002) is shown for comparison. We compare our SN 2008Ddata to the values of the stellar cooling envelope model (dashed 
lines) in Soderberg et al. (2008, S08), using the same parameters (in particular, ^ Rq) well as to the model of Chevalier & Fransson 
(2008, CF08, with their suggested R^, Si9 R©). Right: Bolometric light curves of XRT 080109/SN 2008D (filled circles) computed via two 
different techniques: fitting a BB to the reddening-corrected UVWl Kg fluxes (LbBi top) and via direct integration of the reddening- 
corrected broad-band photometry (Luvwi—'K^, bottom). The black data point at 1.8 day includes the detections in the Far-UV filters 
and thus represents ii7VW2^ifs- We note that the value for the bo lometric magnitude M^^i of SN 2008D will vary depending on the 
passband used over which different authors (808, IMalesani et al.ll2009l) integrate its flux. 



tlie assumption of a BB continuum is expected to break 
down, based on the optical spectra where distinct absorp- 
tion lines develop (see Fig. [5]). While the values of the 
fits are strictly not acceptable except at early times (1 
< At < 5 day), the overall fits are reasonable even be- 
yond this time with no obvious residual trends, and we 
consider the BB measured temperatures, radii, and lu- 
minosities to still be generally quite reliable. Thus, the 
UVOPTIR data are well fit by a BB, especially at early 
times (and At < 5 day), and we do not detect any NIR 
excess (e.g., due to potential dust formation). For clarifi- 
cation, we note that the optical/NIR data to which we fit 
a BB were obtained much later (At > 10^ s) than when 
the main X-ray emission period of XRT 080109 occurred 
{At < 10^ s, § El. For At= 1.8 day, we detect the SN 
in the Swift UVW2 and UVM2 band and include those 
data in the BB fit, though their inclusion docs not have a 
significant effect on the outcome (the parameters change 
by less than 2%). 

We present the resultant BB luminosity evolution in 
Figure [8] (right, top panel) and include th at of SN 1999e x 
for comparison (from Table 8 in Stritzinger "eraI1((200l) . 
for their adopted values of E{B — y) = 0.28 mag, and 
explosion date of JD 2,451,480.5), but with a luminos- 
ity distance recomputed with Hq — 73 km s"^ Mpc"\ 
Both SNe exhibit the initial dip from the cooling stel- 



lar envelope emission, with SN 2008D showing a much 
less pronounced contrast between dip and peak due to 
heating of ^^Ni than in SN 1999ex. 

The second method for computing the bolometric light 
output is by direct integration of the observed broad- 
band photometry. Here we are in the favorable posi- 
tion of having NIR data, as the NIR contribution to the 
bolometric output is usually unknown for SNe Ib/c and 
can be quite large (30-50%) for cer t ain SNe at late time 
(^20_J_avj_ jTominaga et all l2005t iTomita et all 120061 : 
'Modiaz"2007D, a nd even larger for dust-producing SNe 
(SN "2006ic, e .g., lSmith et al.ll2008l: [Pi Carlo et al.ll2008l : 
lModiaz|[200l . 

We integrate the total optical and NIR monochromatic 
fluxes in the region between the effective wavelengths of 
the UVWl and Kg filters via the trapezoid approxima- 
tion. For this purpose, we extrapolate the SED to zero at 
the blue UVWl (2220 A) and red (24,210 A) edges 
of the total set of filters. We plot the resultant bolomet- 
ric light curve (Luvwi^K,) in Figure [8] (right, bottom 
panel) and, for compar ison, also that o f SN 1999ex (based 
onU ^ z integration. IStritzinger et~a l. 2002). 

For SN 2008D, we find that the bolometric luminosity 
based on direct UVWl Kg integration is very simi- 
lar to that based on the BB fits (within < 0.04 dex or 
< 8%), except for the first three data points at At = 
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0.8, 1.8, and 2.8 day (different by a factor 1.32, 1.13, and 
1.09, respectively). SN 2008D reaches bolometric peak at 
At = 19.2 ± 0.3 days with a bolometric peak luminosity 
of logLboi = 42.2 ± 0.1 erg s-\ i.e., A/boi -16.8 ± 0.1 
mag. 

At At = 0.84 day, the peak of the SED is at wave- 
lengths bluer than the observed wavelength coverage and 
explains the missing flux in the direct integration, which 
is the absolute lower limit to the total light output. We 
note that for SN 2008D, the agreement between bolo- 
metric values derived from BB fits and those based on 
direct integration is much better than for SN 1999ex; if 
the emission from SN 1999ex can be as well described 
by a BB as that of SN 2008D, then the true bolometric 
peak luminosity of SN 1999ex is not —17.0 mag (Fig. [51), 
but rather —17.7 mag; i.e., the SN is more luminous than 
previously thoughtm because it was not observed in the 
UV nor in the NIR whose data could have been used to 
construct a full bolometric output. 

In Figure [9l we present the contribution of differ- 
ent passbands to the full bolometric luminosity output 
Luvwi^K, as a function of SN age. At early times, 
the UV is 40% of the bolometric luminosity. At later 
times (30 day after outburst), the NIR is 20%. Most su- 
pernova observations lack these measurements, so their 
derived bolometric flux has significant uncertainties. 
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the temperature of the cooling envelope are weak func- 
tions of the ejecta mass (Mej) and ejecta kinetic en- 
ergy (-Ek), and a stronger function of the radius of the 
progenitor before explosion (i?*) to the power of 1/4. 
Both models assume that the photophere is in the outer 
shock-accelerated part of the SN ejecta density profile, 
but their methods differ. In the following, we use our 
early-time light-curve data to provide constraints on R^, 
if the ratio Ei^/AI^j can be meas ured independently from 
the SN light c urves and spectra (ISoderberg et al.ll2008bl : 
iMazzah et aIll2008HTa"naka et al.ll2009D . 

Analytically, we can rewrite the equations for the ra- 
dius and the temperature of the cooling envelope as given 
in W07 and in CF08 (equations [18] & [19], and equations 
[2] & [5] , respectively) , and generalize such that 
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C = 0.286 for W07 or 3.028 for CF08, 
a = 1.69 for W07 or 1.61 for CF08, 
(3= -0.24 for W07 or -0.27 for CF08, 
where TBB,ev is the BB temperature of the photosphere in 
eV, i?BB = 10^^i?BB,i4 cm, Atd time in days after shock 



breakout, i?K = 10 E'k.si erg, and M; 
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Fig. 9.— Contribution of the UV {UVWl,U: blue circles), 
optical {BVr'i': orange circles), and NIR (JHKs: red circles) 
passbands to the bolometric luminosity Li^t as function of age for 
SN 2008D. The empty circles at 1.8 day denote the contributions 
when including the Far-UV detections (in the UVM2,UVW2 fil- 
ters). 



6.2. Constraining the Progenitor Size of 
XRT 080109/SN 2008D 

Now we compare our data to the fits of lSoderberg et all 

(l2008bD via the coohng BB envelope model of 

Waxman et al .l (|2007|), an d to the predictions of 



Chevalier & F ranssonr(|2008f ) as the source for the early- 
time emission, where the stellar envelope expands and 
cools after the passage and breakout of the shock wave. 
In both models, the time evolution of the radius and 



Thus, the derived BB fits to the early-time data yield 
independent measurements of Ri, for each epoch, for a 
given Ek/Mc}, whose uncertainty does not influence the 
value of significantly, given the very weak dependence 
of R^, on Ex/Mcy While the power-law scalings of the 
parameters are fairly similar for the two models, their 
largest discrepancy lies in a factor of ~10 difference in 
the constant C, which is primarily caused by the factor 
of ^2 difference in the temperature equations, as noted 
by CF08, and by the high sensitivity of i?* on Tbb. 

Here we explicitly calculate R^, for the two models 
with our own data combined with published values of 
- Ei^/Mej and include all relevant sources of uncertain- 
ties; using the slightly different values for Ek/Mb] for 
SN 2008D fro m the literature (ISoderberg et al.l l2008bl: 
IMazzah et all 120081: iTanaka et al.ll2009D . and our esti- 
mates for Tbb and i?BB (Table |8|), we compute i?* for 
lOO.Oeach of the first four epochs (0.8 < At < 3.8 day), and for 
the range of adopted extinction values {E{B — F)Host — 
0.6 ± 0.1 mag). Taking an average over the epochs and 
using the different reddening values as bounds, we obtain 
for the W07 model i?W07 ^ (8.6±crEx + crE/M) x 10^" cm, 
where ctex and CTe/m refer to the systematic error due 
to uncertain extinction and uncertain E^/Mcj, respec- 
tively, with ctex — 4.8 and cte/m = 0-4. In other words, 
fjWor ^ i.2±O.7±O.O6R0. Likewise for the CF08 model, 
i?CF08 ^ (86±45±4) X 10^° cm = 12±7±0.6 Rq. Thus, 
Ri, as derived from the cooling stellar envelope emission 
is 10 times larger according to the CF08 model than ac- 
cording to the W07 model. CF08 had noted this fact, and 
furthermore, that their larger progenitor size is consistent 
with R^~'^'^^ « 9 R0 , based on their simple interpreta- 
tion of the published X-ray spectrum of XRT 080109 as 
a purely thermal X-ray breakout spectrum. 

In addition to the literal application of both models, 
it is important to keep in mind the caveats and simpli- 
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fications that enter them, as both models are consistent 
with the basic idea of emission from a post-breakout cool- 
ing stellar envelope. Any simplifications that lead to a 
slight mis-estimation of Tbb,ov will lead to a large effect 
onto Ri, due to its large temperature-sensitivity. Caveats 
here include equating effective temperature (as given by 
the models) with observed color temperature, and the 
fact that for temperatures during the time of observa- 
tions (10"* K), the opacity is probably not set by simple 
Thomson scattering as assumed in the models, but is 
composition dependant (E. Waxman, 2009, private com- 
munication). These aspects mark areas of improvement 
that need to be addressed if measurements like the ones 
for SN 2008D are to be used to routinely and precisely 
measure SN progenitor sizes in the future. 

For illustration, we overplot in Figure [8] the predicted 
values f or the Waxman model us ing the parameters put 
forth bv lSoderberg et al.l (|2008bD for SN 2008D, namely 
= 10" cm, Mej = 5 Mq, and = 2 X lO^^ erg 
(see their caption to Figure 3). Our values are consistent 
with their fits, within the uncertainties. For compari- 
son we also plot the predic ted values for the model of 
IChevalier fc FranssonI (j2008l ) with i?^ = 9 Rg . 

We can now directly compare our estimated i?^ 
of XRT 080109/SN 2008D to stellar radii of Wolf- 
Rayet (WR) stars, which posse s s den se winds (see 
lAbbott fc Cmtil [19871 : iCrowthed [20071 for reviews), 
and are the mos t likely progen itors of stripped- 
envelope SNe (e.g., lFilii)Denkolll991h, includin g that of 
XRT 080109/SN 2008D (iSoderberg et al.ll2008bl ^. Specif- 
ically, bona fide WN stars — the nitrogen-sequence spec- 
tral subtype of W R stars with their heli um layer intact 
but without H (see lSmith fc Contil[2008l for the distinc- 
tion between WN and WNH, i.e., H-rich WN stars) — 
are the likely progenitors of some SNe lb. Their radii are 
in t he range 1-20 R^ with ty pical sizes between 5—10 
R,7. jHerald et al.l [200lt iHamann et all 120061 : ICrowthed 
I2007D . Provided that those models and their data are 
accurate, w 0.5-1.9 Rq as suggested by WO 7 is 

smaller by a factor of 5—10 than typical values for such 
WN radii, while the size of R'^^"^ = 5 - 20 Rq as sug- 
gested by CF08 is more in line with data of typical WN 
stars. 

Independent of t he considerations above, 
ISoderberg et al.l (|2008b[ ) infer from their extensive 
radio observations that the progenitor had a steady 
wind with a mass-loss rate consistent with those of WR 
stars. 

7. COMPARISON WITH OTHER SUPERNOVAE 

The optical signature of postshock breakout, namely 
a cooling stellar envelope (extremely blue color and 
dip in the first few days after explosion) has been ob- 
served in two stripped-envelope normal SNe: SN lib 
1993 J ([Schmid t et al." 19931 [Richmond et~aI1 Il994f ) and 



SN lb 1999ex (Stritzin ger et al.l 



200^. Other SNe with 



early-time observations (at At w 1 day or earlier) in- 
clude the broad-lined SN Ic 2 006a.i connected with low- 
luminosity GRB/XRF 060218 ([Campana et al.l[2006L but 
see § m, GRB 980425 /SN 1998bw (for the V and R 
filters.TGalama et all [1998[ ). the famous SN II 1987A 
([Hamuv et al.[ [19881 ). and most recently, two SN II 
found by the Supernova Legac y Survey and de tected 
with the UV satellite GALEX '([Schawinski et al.[ |2008[ : 



iGezari et al1[200l . 

7.1. Photometry: Cooling Stellar Envelope and 
Radioactive Decay as Power Sources 

Figure (TOl compares the optical light curves of 
SN 2008D to those of other SNe with either observed 
shock breakout or subsequent phase of cooling stellar en- 
velope emission, where the SNe are offset to match the 
peak magnitude of SN 2008D in the respective bands. 
Careful comparison with other stripped-envelope SNe 
caught shortly (~ 1 day) after shock breakout (SN 1993J, 
SN 1999ex) reveals that while the light curves (and spec- 
tra) of SN 2008D and SN 1999ex are similar during the 
^^Ni-decay phase, the light curves related to the cool- 
ing stellar envelope are different among these three SNe. 
Furthermore, we note that the B — V color evolution 
of SN 2008D is very similar to that of SN 1999ex dur- 
ing the time of ^^Ni decay, but very different during the 
shock-breakout phase (Figure [TT|l . The color curve of 
SN 1999ex has a steep rise to and decline from the red 
during 1< At < 10 days, peaking at B — V=l mag, while 
that of SN 2008D evolves gradually during this time and 
reaches B — V=0.5 mag at most. 

We leave it as a future topic to model the diversity 
of the very early-time light curves of SNe 1993J, 1999ex, 
2006aj, and 2008D in terms of progenitor radius, £'k, and 
Mcj , and to compare these values with those derived from 
other methods (e.g., light curve and nebular models). 
We note, however, that the temporal evolution of the 
inferred BB temperature and radius of SN 2008D are 
quite diffe rent from those d erived for SN 1999ex (Fig. [H 
left panel. [Stritzinger et al.[ r2002). 

Comparing the light-curve shapes of the different SNe 
with respect to date of V^-band maximum, we find that 
SN 2008D declines more slowly than any of the other SNe 
and has a wider peak. One route of parameterizing the 
SN hght-curve shape is using the Ar?ii5( X) parameter, a 
formalism initially developed for SNe la l|Phillipsl[1993f l. 
which is defined as the decline (in magnitudes) during the 
first 15 days after maximum brightness in the passband 
X. For SN 2008D, we measure AmisiU) = 1.2 mag, 
Ami^{B) — 0.8 mag, and Ami^{V) = 0.6 mag. 

7.2. Optical Spectroscopy: Emergence of He I Lines 

The very early spectra [At « 3 — 5 day) actually re- 
semble those of S N 2006a.i ( sec top two spectra in Fig- 
uredH see also iMazzah et alj, 2008: Malcsani et al. 2009), 
the SN connected with GRB/XRF 060218, which ini- 
tially led to the ide ntification of this SN as a b road- 
line d SN ([Blondin et al. 20081 [Kong et al.|[2008aD (but 
see iMalesani et al.l [2008l) . This spectral resemblance 
suggests that in both systems, the material constitut- 
ing the photosphere at early times is at high veloc- 
ity ([Mazzali et al.[ [20081 ). as expected from homologous 
expansion. Indeed, the velocity we determined via 
SYNOW, Uphot « 30,000 km s"^ at At = 1.8 day in 
the spectrum with the double-absorption feature (see 
§5.3p . is similar to those of broad- lined SNe Ic at the 
same epoch measur ed in a slightly different manner 
([Mazzali et al.[[2008[ ). 

However, while SN 2006aj and other broad-lined SNe Ic 
continue to have a significant amount of mass at high ve- 
locity (at least 0.1 Mq at 20,000-30,000 km s"^ of the 
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Fig. 10. — Comparison of the optical Swift/VYOT (filled circles) and ground-based (empty circles and triangles) light curves of SN 2008D 
with those of other normal SNe (left) and SNe associated wi th low-luminosity GRBs {right, but sec § [9ll, for which shoc k breakout was 
obs erved or the time of explosion is well known. Data are from'R ichmond e t al.! fl994 } for SN1993J , .Strit zinger et al.l (12002 ) for SN 1999ex, 
[Galama et al. (1998) for SN 1998bw, and Camoana et al. (2006) (UBV) and Modi az et al.l (|200g i {r'i') for SN 2006aj. The offsets refer 
to the UBVRt light curves of SN 2008D, while its r'i' light curves are offset by r' - 2.6 and i' - 4.3 to match RI. The rest of the SNe 
are offset to match the peak magnitude of SN 2008D in the respective bands. For SN 199 8bw and SN 2006aj, the SN explosion date 
was taken to be coincident with the onset of GRB 980425 and XRF 060218, respectively (Gal ama et al.| [r998l; ICampana et al.ll200(3l '). For 
SN 1993J and SN 1999ex, the explosion dates were estimated from optical nondetections on the night before (for SN 1999ex, JD 2,451,480.5; 
IStritzing er ct al. 2002'). and on the night during discovery (for SN 1993J, JD 2,449,074.5 ± 0.05; Wheeler et al. 1991). 

In SN 2008D, the He I lines become apparent as early 
as Ai « 6 day (i.e., iymax ~ — 13 day) and are fully 
identified at At « 11.7 day (tymax ~ —6.9 day). Fig- 
ure [T3] reveals that the velocities span a large range for 
normal SNe lb; e.g., at maximum light, they range from 

-14,000 km s~^ (for SN 2004gq) to 10,000 km s"^ 

(for SN 2005hg), and seem to fo llow a power-law decline. 
For comparison, [Branch et al.l ([2002 ) showed that the 
He I velocities in the SNe lb of their limited sample fol- 
lowed a standard pattern and traced fairly closely (within 
2000 km s^^ for the same epoch) the photospheric veloc- 
ities derived from synthetic fits to Fe II lines, which we 
overplot in Figure [T3] (dashed line, using their power- 
law fit fphot t^2/(n-i) — 3.6). However, the 
SNe lb in our sample show a larger departure (up to 
6000 km s~ ) from the corresponding standard photo- 



total of 2 M(T) ejec ted mass as derived from models in 
iMazzali et al.l 12006). at At « 12 day, the receding pho- 
tosphere in SN 2008D reveals the lower-velocity layers 
and the development of He I lines which SN 2006aj and 
other broad- lined SNe Ic do not show (Fig. fT^ . Spectral 
synthesis modeling yi elds an estimate of ^0.03 Mfr^ at 
y > .1c for SN 2008D f Mazzah et al. ''2008:'Ta Iiaka et all 
[2009f) . a factor of 10-50 less than for GRB-SNe. 

As Figure [T^ shows, the later spectral evolution of 
SN 2008D resembles that of a normal SN lb, where the 
He I lines are as pronounced as in SN 2005hg and more 
so than in SN 1999ex. 

We now investigate the temporal behavior of the He I 
lines, which are the spectral hallmarks of SNe lb. In 
SN 2008D and in a larger sample of SNe lb, we deter- 
mine the blueshifts in the maximum absorption in the 
He I lines via the robust techniq ue of meas uring SN line 
profiles as presented by Blondi n et all (|2006) . 

We refer to our velocity measurements as "He I 
absorption-line velocities" in the rest of this work. In 
Figure [13] (see also Table [9]) we present the measured 
values for He I A5876, the strongest optical He I line in 
SN 2008D, along with ab sorption-line v elocities of a sam- 
ple of other SNe lb fsee lModiazll2007[ ). measured in the 
same fashion. 



spheric veloci t y, larg er than those of SN 1998dt, which 
[Branch et al.l (|2002D had declared as the exception in 
their sample. In addition, both SNe 2008D and 2005hg 
show He I line velocities that are lower than the cor- 
responding photospheric velocity for dates before maxi- 
mum light. This might mean that the line-emitting ma- 
terial is located at radii smaller than the bulk of the 
material that makes up the photosphere, or that there 
are optical-depth effects at play. 
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Fig. 11.— Comparison of the S - V color curve of SN 2008D 
dereddened by E{B ~V) = .6 mag, with those of o ther SNe. The 
B-V color of SN lib 1993J IIRichmond et al.|[T99i) is dereddened 
by E{B — V) = 0.08 mag (we adopted the lower value derived by 
Richmond et al. and not their higher alternative of E{B — V) = 
0.28 mag), and t hat of SN lb 1999cx by E{B - V) = 0.3 mag 
llStritzineer et al .' 2002). 

It has been long suggested that the He I lines are due to 
nonthermal excitation and ionization of gamma-ray pho- 
tons produced during the radioactive decay of ^^Ni and 
56Co (IHarkness et alJUOSTl: ILucvI119911 IMazzah fc Lucvl 
I1998D . The fact that the gamma-ray photons need time 
to diffuse through the ejecta to reach the He layer has 
been invoked to explain the late emergence (after max- 
imum light) of He I in some SNe lb (e.g., iFilippenkol 
Il997l and references therein). A challenge for those who 
study the complicated nonthermal excitation effects in 
SN line formatio n (Baron et al. 1999; Kasen et al. 200j; 
iDessart fc Hillieill2008i) would be to self-consistently ex- 
plain for SN 2008D its early emergence of He I lines (~5- 
6 day after explosion, and ^12-13 day before maximum 
light) by excitation of the relatively small amount of 
s ynthesized ^^Ni (0 . 05-0.1 M^t^ , ^S oderb erg et al...2008bl : 
[Mazzah et all 120081: iTanaka et alj|2009f) . 

7.3. NIR Spectroscopy as a Possible Discriminant 

Figure [Mj shows the four NIR spectra of SN 2008D 
(taken with the IRTF) in which we detect five prominent 
blueshifted He I lines. We tabulate their velocities as 
inferred from maximum absorption in Table O 

Comparison of the NIR spectra of proper SNe lb 
1999ex, 2001B, and 2008D with those of SNe Ic (see 
iTaubenberger et al.l 12006 and references therein) shows 
that all SNe display a strong feature at ^^1.04 but 
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Fig. 12.— Comparison of the MMT and APO spectra of 
SN 2008D taken at At = 3.74 day and 11.66 day, respectively, 
with those of SN 2006aj at similar ages (At = 4.95 and 12.85 day, 
respectively), and with those of other SN lb and SN Ic. SN 2006aj 
was the SN associated with XRF 060218 l IModiaz et al.ll2006l) . 
The early-time spectra of other SNe lb and I c include, dered- 
dencd by the indicated amounts: SN 1999cx fStrit zinger et aLl 
|2002; Hamuv ct al. 2002), SN 2005hg (Modjaz 2007), and SN Ic 
19941 (FilipDcnk o et al.lll995IV The optical scries of He I lines, 
defining SN 2008D as a SN lb, are indicated (namely. He I AA4471, 
5876, 6678, 7061). We note the similarity between SN 2008D and 
SNe 1999ex and 2005hg. 
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Fig. 13. — Temporal velocity evolution of He I A5876 in SNc lb: 
SN 1998dt ( Mathcson ct al. 2001), SN 2005bf (Tominaga ct aT] 
[20051 : lModiaa,2007) . SN 2004gq and SN 2005hg ( Modjaz 2007|)7 
and SN 2008D (this work). He I A5876 is usually the strongest 
He I line. The alssorption velocities decrease over time, except for 
the maverick SN lb 2005bf (see text). The velocities for SN 2008D 
lie within the range of observed values, with the earliest measure- 
ments before maximum for a normal SN lb. The thick dashed line 
is the functional form for the photospheric velocity as fit by Branch 
et al. (2002, B02) to SN lb data. The thin dashed lines demarcate 
the extent in their uncertainty. See text for details. 

only SNe lb show additional absorption at 1.95—2.0 /im 
due to blueshifted He I 2.058 /xm, while none of the 
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Fig. 14.— Three NIR spectra of SN 2008D taken with the IRTF. 
The dates of observation are indicated, along with SN ages (At) 
with respect to time of shock breakout (in days). The five detected 
He I lines are indicated by their rest wavelengths. The spectra are 
shown as both raw and Fourier smoothed, and are normalized and 
offset by arbitrary amounts for clarity. Note the strong presence 
of blueshifted He I 1.0830 at ~ 1.02-1.05 /xm and blueshifted He I 
2.0581 fj,m at ~ 1.97—1.98 /im. The broad absorption at 1.05 ^m 
probably includes an additional blend of other lines (mostly Fe II, 
but also Mg II, C I, and Si I). 

SNe I c do. For the few SNe lb in general (jGerardv et al.l 
12004 ) and for SN 2008D in particular, the absorption 
velocity of the strong blueshifted He I 2.058 /xm is very 
similar to the strong blueshifted He 1 1.083 /xm if we iden- 
tify the redder absorption trough of the multi-component 
feature at 1 fim as He I 1.083 /im. As Figure [T5l shows, 
the blueshifted He I 2.058 fim line is narrower and more 
well-defined than the broad blend at 1 /im. 

In conclusion, SN 2008D gives supporting evidence to 
the argument that the broad absorption at ^1.04 fim 
is a poor diagnostic for helium, since it is most likely 
a blend of F e II, C I, Ca II, a nd S I at the same 
wavelengths (iMillard et al.l 119991 : iGerardv et all 120041 : 
iSauer etahl 120061 ) and He I in SNe lb. However, claims 
of significant He in optically classified SN Ic rely on iden- 
tifying the same strong absorption featu re at ^1.04 /im 
exclusively with He I 1 .083 /zm ( Fihppe nko et al.l 119951 : 
IClocchiatt i et al.l ll996': Pat at et al.l l2001). Our study of 
SN 2008D implies that determining unambiguously if in- 
deed helium is present in SNe normally classified as Type 
Ic requires a iiT-band spectrum to test for the presence 
of the clean and unblended He I 2.058 /im feature. 

8. DOUBLE-PEAKED OXYGEN LINES AND ASPHERICAL 
EXPLOSION GEOMETRY 

In Figure [TBI we present the Keck spectra of SN 2008D 
taken at tvmax = 91.4 day (i.e., At = 109 day) on 
April 28, 2008, and iymax = 132.4 day, on June 7, 2008. 
They show nebular emission lines on top of some resid- 
ual, underlying photospheric spectrum, which is more 
pronounced in the earlier spectrum. This is expected, 
as SNe appear to turn fully transparent ~200 day after 
maximum light, and some SNe take up to 1 year to com- 
plete the transition (Mazzali et al. 2004). Here, we are 
witnessing this transition period for SN 2008D. Below we 
analyze the nebular emission lines of SN 2008D. 

At sufficiently late times, the SN ejecta become fully 
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Fig. 15. — The two main helium lines detected in NIR spectra 
of SN 2008D in velocity space for two epochs: 2008 Jan. 21.51 
(At = 11.95 day; solid Unes) and Jan. 28.51 (At = 18.95 day; 
dashed lines). The top plot shows the NIR spectra with respect to 
He I 1.083 /im as zero velocity, and the bottom one with respect 
to He I 2.058 /im. The largest difference between the two lines is 
the effect of blending of other elements with the He I 1.083 /im 
line. The blueshifted He I 2.058 fira line is narrower and more 
well-defined than the blend around 1 fira, of which we take the red 
absorption trough to be due to He I 1.083 /tm. Both He I lines have 
similar velocities based on the blueshifted maximum absorption 
(see Table |9]l. The spectra are normalized to the local continuum. 

optically thin in the continuum; hence, spectra obtained 
during this time period afford a deeper view into the 
core of the explosion than spectra taken during the early 
photospheric phase. In the absence of hydrogen, oxy- 
gen is the primary coolant in the ejecta of stripped- 
envelope SNe at late epoch s when the gas is neutral 
or at most sing ly io nized (jUomoto fc Kirshneil Il986l : 
iFransson fc Chevaheilll987l ). and when densities are suf- 
ficiently low for forbidden lines to be the strongest ones 
in the spectrum. In SN 2008D we detect [O I] AA6300, 
6363 and O I A7774 along with [Ca II] AA7291, 7324 and 
the strong Ca II NIR triplet. 

To study the line profiles of SN 2008D more closely, 
we plot the main relatively unblended emission lines 
([O I] AA6300, 6363, O I A7774, and [Ca II] AA7291, 
7324), in velocity space in Figure [iTl The two oxygen 
lines show the same unusual double-peaked profile, while 
[Ca II] AA7291, 7324 does not exhibit it. For both oxygen 
lines, the redder peak is at zero velocity, and the trough 
between the two peaks is shifted by ~ —800 km s~^ 
from zero velocity. The total separation between the 
two peaks is 2000 km for both the forbidden and 
the permitted oxygen lines. With the chosen zeropoint 
for [Ca II] AA7291, 7324 (namely 7307.5 A, the straight 
mean of the two calcium lines) , its peak is also blueshifted 
by ^-^800 km s~^, but blending with other lines ([O II] 
AA7319, 7330, Fe) could be affecting the line shape and 
the exact determination of zero velocity. 
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Fig. 16. — Spectra of SN 2008D at 3 and 4 months after maximum 
light (April 28, and June 7, 2008). The main emission lines in the 
spectrum are marked, and their rest wavelengths indicated. We 
caution that the line at 7300 A could be a blend of [O II] AA7319, 
7330 and some Fe lines, in addition to [Ca II] AA7291, 7324, as 
marked. 

While optical-depth effects sometimes generate double- 
peaked [O I] AA6300, 6363 because of its doublet na- 
ture (most no tably in SN 1987A: ISpyromiliol 119911 : 
lLeibundgut"era l. 1991; iLi fc McCravl I1992D . it is not 
the culprit in SN 2008D, since (a) the separation of the 
doublet (3000 km s~^) is larger than the observed sep- 
aration between the double peaks (2000 km s~^), and 
(b) O I A7774, a multiplet with very small separations 
(at most 3 A, i.e., 130 km s~ , smaller than our spec- 
tral resolution), shows the same double-peaked profile as 
[O I] AA6300, 6363 in SN 2008D. The [Ca II] AA7291, 
7324 line does not show the same double-peaked profile, 
but this may well be because the emission comes from 
smoothly distributed, preexisting (nata l) calcium simi- 
lar to what was seen in SN lib 1993J (|Matheson et al.l 
I2000al) . 

Such double-peaked oxygen emission lines have 
recently been ob served in a n umbe r of stripped- 
envel op e SNe (llVlazzali et al.l 120051: llVIaeda et all 

20071: iModiaz et al.l l2008cl: iMaeda et al] [200a 

Taubenberger et"aI1 l2009 f). and have been interpreted 
in the framework of global asphericities, where the 
emitting oxygen-rich ejecta are situated in a ring- or 
torus- hke structure v iewed along the equatoria l plane . 
Both llVlodiaz et all (|2008cD and iMaeda et all (|2008D 
suggest that such global asphericities are common 
during core collapse. Macda et al. (2008) specifically 
suggest that their observed fraction of double-peaked 
profiles (40 ± 10%) is consistent with the hypothesis that 
all stripped core-collaps e events are mildly aspherical, 
given their models in IMaeda et al.l ()2006[ ). whereas 
iTaubenberger et all (120091) find a lower observed frac- 
tion (5-18%). Recently. iMilisavlievic et al.l (|2009D have 
argued that it is difficult to explain all line profiles 
observed in their recent SN sample via purely geometric 
effects. However, their conclusions only concern SNe 
whose double-peaked profile has a velocity separation 
close to that of the doublet [O I] AA6300, 6363 and 
does not show it in lines other than the doublet; those 
features are not observed in SN 2008D (see Fig. [16] and 
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Fig. 17. — The main late-time emission lines of SN 2008D in 
velocity space in the rest frame of NGC 2770. The zeropoint for 
[O I] AA6300, 6363 here is 6300 A while that of [Ca II] AA7291, 
7324 is taken at 7307.5 A (straight mean of the two calcium lines). 
We caution that [Ca II] might be a blend with [O II] and Fe. 
Zero velocity is indicated by the dashed line. The two oxygen 
lines show the same conspicuous double-peaked profile, indicat- 
ing that the two horns cannot be due to the doublet nature of 
[O I] AA6300, 6363. In both cases, the redder peak is at zero ve- 
locity, and the trough between the two peaks is shifted by ~ —800 
km s~^ from zero velocity for both spectra. With the chosen zero- 
point for [Ca II] AA7291, 7324, its peak is also blueshifted by ~800 
km s~^, though blending could be affecting the line shape and the 
exact location of zero velocity. See text for more details. 

above) . 

While other SNe lib and lb have been shown to ex- 
hibit some small-scale structure in oxygen, indicating 
some clumping in the den sity distribution of oxygen 
probably due to instabilities (Filiopenko & Sargent"1986|: 
ppvr omilio 1994; SoUerman et al. 1998; IMathcson ct al] 
l2Q0Qal: lElmhamdi et al.ll2004D . the clear symmetric dou- 
ble peaks in SN 2008D are different and may rather 
indicate a global scale for the anisotropy. The mech- 
anism giving rise to such a global anisotropy could 
lie in the intrinsic explosion physics of core collapse, 
since each of the various t heoretical core-coll apse mod- 
els, be it neutrino d riven (jScheck et al.ll2006l). acoustic 
(iBurrows et al.ll2006r). or magneto driven ((Burrows et al.l 
l2007HDessart et al.ll2008 l) , finds that symmetry breaking 
is essential for a successful explosion. 

Alternatively, in case the progenitor of SN 2008D is 
part of a binary, binary interaction or a merger might 
be modulating the geometry of the explosion, though 
this speculation will be very hard to verify. Neverthe- 
less, other lines of evidence, suc h as polarization stud- 
ies o f SNe II, lb, and Ic (e.g., iLeoi iar d fc Filippcnkol 
[2005I) . neutron-star kick velocities (jWang et al. ,2006i ), 
and young SN remnant morphologies (e.g., Cas A; 
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iFesen et al.ll2006l ) suggest that asphericities are generic 
to the core-collapse process. 

We note the blueshift from zero velocity (of the order 
~ 800 km s^ , Figm-c [T7)l in the oxygen trough, if the 
trough between the two peaks is regarded as the sym- 
metry center. The same blueshift might be also present 
in the calcium peak, depending on the choice of zero 
velocity, which is more uncertain than for the oxygen 
lines due to possible blending with other lines. This 
blueshift, if real, could be caused by residual opacity be- 
fore r eaching the full nebular phase (jTaubenberger et al.l 
l2009l ). pecuhar distribution of ^^Co which is powering 
the emission at such later epochs, large-scale asymme- 
try, or offset of the SN rest-frame velocity from that of 
NGC 2770 by 800 km s~\ However, given that the ro- 
tation curve of th e host galaxy, NGC 2770, extends to 
only ~400 km s~^ ()Havnes et al.lll997[ ). the last option is 
least likely. 

Future multi-epoch late-time spectra (especially those 
taken at tvmax ^ 200 day) are needed to differentiate 
between the various scenarios as each one predicts a dif- 
ferent behavior of the blueshift over time. If it is due 
to residual opacity, the most probable case, we predict 
that the blueshift will go to zero for spectra at iymax ^ 
200 day. We can exclude dust formation as the culprit, 
since the apparent blueshift does not increase over time, 
contrary to the behavior of bl ueshifted lines in S N where 
dust had been formed (e.g., ISmith et al1l2008f ). If the 
blueshift con tinues to remain const ant over time (as seen 
in SN 1996N: ISollerman et al.lll998[) . then bulk blueshifts 
or unipolar blob ejection s might need to be invoked as 
suggested for SN 2005bf (|Maeda et al.ll2007l ). 

In summary, we infer from the presence of double peaks 
in the two different oxygen lines in SN 2008D that global 
asphericities were probably involved during the explosive 
event of SN 2008D, regardless of the specific mechanism 
that actually leads to a successful explosion. 

9. X-RAY INTERPRETATION 

emissi on of XRT 80109 
(l2008bD and 



The origin of the X-ra' 
IS highly debated. ISod erbcrg et al.l 



I Chevalier fc FranssonI (|2008i) attribute it fully to the pro- 
cesses relate d to core- collapse s hock-breakout emission , 
while others (|Xu et al.ll2 008: Li 2008 HMazzah et al.l2008f l 
claim it as emission caused by a stifled or weak GRB jet. 
[Sod crbcrg ct al. (2008b) adopt the pure p ower-law X- 
ray sp ectrum fit and use the calculations of I Wang et al.l 
()2007D to explain the power-law shape in terms of upscat- 
tering (i.e., Comptonization) by the circumstellar wind 
medium of thermal shock breakout photons with cner- 
gies < 0.1 keV. On the other hand. Chevalier & Fransso3 
12008') favor the pure blackbody fit and explain the arti- 
ficially small BB radius as the product of a scattering at- 
mosphere, where the emission is generated in a relatively 
deep layer but effectively emitted at a larger radius. Even 
for XRF 060218/SN 2006aj, which had a much higher 
S/N spectrur n, the exact origin of the X-ray emission is 
stiU debated (ICameana et alJ2006l:fGhisellini et al.ll2007l : 
ISoderberg et al.l 120061 : 1 Waxman et alll2007| ). 

While the above authors attempt to test the sys- 
tem of XRT 080109 for emission from relativistic ejecta, 
their conclusions differ and are heavily model depen- 
dent. In the next section, we explore a phenomenolog- 
ical approach where we compare the observed proper- 



ties of XRT 08010 9 with those of so-called X-ray flashes 
(jHeise et al.ll200i[ l. high-energy transie nt that are likely 
softer analogues of cosmological GRBs (iSakam oto et al.l 
I2008f ). We note that recently, iBietenholz eTa l. (200^ 
find no evidence for a relativistic, off-axis jet in SN 2008D 
based on their VLBI obs ervations of SN 2008D. Very re- 
cently, iKatz et all ()2009[ ) have argued that non-thermal 
X-ray spectra may be common during the breakout of 
radiation-mediated non-relativistic shocks, as suggested 
to be the case for XRT 080109. 

9.1. Comparison with X-ray Flashes 

Here we compare three properties of XRT 080109 with 
those observed in a sample of classical XRFs: duration 
(as parameterized by Tgo), the shape of the X-ray spec- 
trum, and X-ray-to-optical flux ratio. In addition we 
discuss its energy output. 

While in each case XRT 080109 is not necessar- 
ily an obvious outlier in the observed distribution, to- 
gether all three considerations paint a strong picture that 
XRT 080109, along with XRF 060218/SN 2006aj, might 
be different phenomenologically from the rest of classi- 
cal XRFs. However, future studies trying to uncover the 
inherent distribution will need to carefully consider de- 
tector effects and rest-frame properties. 

1) Duration: XRFs observed by HETE-2 and Swift are 
seen to have similar temporal profiles and durations com- 
pared to long-duration GRBs, exhibiting a FRED (i.e., 
fast rise exponential decline) profile. XRT 080109 can be 
adequately fit by such a FRED profile (§ 19. 1|) . but the 
transient's duration is an order of magnitude longer than 
the median Tgo for other XRFs and X-ray-rich GRBs 
(XRR s), which is roug hly 30 s ([Sakamoto et all I2005L 
l20f)8h . Both XRT 080109 and XRF 060218/SN 2006aj 
have the longest durations (Tgo = 470 ± 30 s and 2100 
±100 s, respectively). We note the caveat that satellites 
are typically inhibited from triggering on long timescales 
due to background modulations, and XRT 080109 could 
not have been detected with any other survey X-ray tele- 
scope before Swift. 

2) X-ray spectrum: XRF 060218/SN 2006aj has a 
significant X-ray component that (a) is thermal and 

b) appears to soften over t i me, w ith very high S/N 
Campana et "all 120061 : lButleill2007ar ). The X -ray spec- 
trum of XRT 080109 is of much lower S/N, and is 
well-fit by a combined PL-I-BB (§ [2]). While the X- 
ray emission of XRF 060218 fades at a rate compara- 
ble to that often observed for GRB X-ray afterglows, 
its sp ectrum is markedly softer (e.g., ISoderberg et all 
120061 ) and also softens strongly in time (lButledl2007'air 



This is in contrast to GRB X-ray afterglows, which ex- 
hibit little spectral evolution at the same epoch: GRB 
and XRF X-ray afterglows have photon indices of F w 
2, while th at of XRF 060218/SN 20 06aj was F = 
3.41 ± 0.13 (|Butler fc Kocevskil 120070 ). There is evi- 
dence that the XRT 080109 late-time X-ray spectrum 
(from the Chandra data) was also softer with F = 3.6l^'4 
(|Poolev fc Soderberg|[2008l : ISoderberg erani2008bf ) than 
the early-time spectrum. 

3) We also compare the radiation output in 
XRT 080109 for different wavelengths to that of XRFs 
and their afterglows. 

At the dates of our optical and NIR photometry (at 
t-to^ At = 0.84 and 0.71 day, i.e., 70 and 60 ks), the 
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X-ray luminosity is in the range ^ 2 x 10'^^ '^^ erg s \ 
Thus, the X-ray luminosity (for 0.3-10.0 keV) is 3-4 
orders of magnitude lower than the luminosity inferred 
from the optical observations (see ij lS.ip . thus, ix/^opt ~ 
10-3 - 10-4 at 70 ks for XRT 080109/SN 2008D. This 
corresponds to a ratio of fluxes /o/x = 4 x (10^-10'^), 
while the observed ratio is /o/x < 70 for classical XRF 
aftergl ows (D'Alessio e t al.l |2006[) and GRB afterglows 
(|De Pasa uale et al. 200^ at comparable times (using the 
same definition and units of /o/x as in iD'Alessio et al.] 
|2006[) . In other words, the optical to X-ray ratio for 
XRT 080109/SN 20G8D is a factor of - 10^ larger than 
for XRF afterglows. For XRF 060218/SN 2006aj, the 
ratio cannot be precisely determined since it was not ob- 
served in the desired bandpass at such early time. 

The two main features commonly observed in XRFs 
are not measured for XRT 080109, namely detection of 
gamma-ray emission and peak of the energy distribution 

(-Epoak)- 

In summary, XRT 080109, along with 
XRF 060218/SN 2006aj, appears qualitatively dif- 
ferent from the observed distribution of classical XRFs, 
when considering all three discussed properties together 
(duration, spectral evolution, and optical-to-X-ray ra- 
tio). Finally, the prompt high-energy emission released 
in XRT 080109 iEx,iso « 6 x 1 0^^ erg) is 3-4 orders o f 
magnitude lower than in XRFs (jSoderberg et al.ir2008bf ) . 

10. CONCLUSIONS 

We have presented extensive early-time photometric 
{UVW2,UVM2,UVWl,UBVRr'Ii'JHK,) and spec- 
troscopic (optical and NIR) data on SN 2008D as well 
as X-ray data analysis on the associated XRT 080109. 
Our data span a time range of 17 hr to 35 day after 
the detection of the X-ray transient and establish that 
SN 2008D is a spectroscopically normal SN lb at maxi- 
mum light, with a relatively long rise time (18 day) and 
a modest optical peak luminosity of My = —17.0 ± 0.3 
mag and Mb = —16.3 ± 0.4 mag, after correcting for 
EiB - F)Host = 0.6 ± 0.1 mag. 

We also present the earliest NIR (JHKs) detections 
(at 0.71 day). Our densely time-sampled optical spec- 
tral sequence, along with our NIR spectra, uncovers the 
emergence epoch of the helium lines. The early-time 
spectra reveal high expansion velocities (15,0000—30,000 
km s~^) up to 6-7 day after outburst, and one specific 
spectrum (at At — 1.84 day) shows a transient, promi- 
nent double-absorption feature at blue wavelengths (~ 
4000 A) that is no longer visible in spectra taking 1 
day later. In spectra taken 3 and 4 months after maxi- 
mum light, we detect double-peaked profiles in the oxy- 
gen lines that we take to indicate global asphericities in 
the SN core ejecta, and argue that they cannot be due 
to optical-depth effects. Furthermore, we recommend for 
future searches for helium in SNe Ic that J- and iiT-band 
spectra are obtained to test for the presence of both He I 
1.083 ^i^Jl and He I 2.058 /xm. 

The comprehensive data presented here allow us to 
construct a reliable measurement of the bolometric out- 
put for this stripped-envelope SN, including UV and NIR 
contributions that are usually lacking for other SNe. For 
data At > 3 days, we find that the bolometric luminosity 
based on direct UVWl Ks integration is very similar 



to that based on the blackbody fits to the broad-band 
photometry (within <0.04 dex or < 8%), such that the 
UV-NIR emission of SN 2008D may well be described 
by that of a black-body. We furthermore compare ou r 
early-time light curves to fits by Soderberg et al. (2008b'), 
who use the calculations of Waxman et al . (2007) for the 
cooling stellar envelope blackbody model, and test the 
predictions of the alte rnative cooling envelope model of 
IChevalier fc FranssonI |2008). Combining our compre- 
hensive data with estimates of £'k and A/cj from the 
literature, we estimate the stellar radius Ri, of its prob- 
ab le Wolf-Rayet progeni tor. Acc ording to the models 
of IWaxman et al.f (|2007l ) and of IChevalier fc FranssonI 
(l2008l) . we derive i?f °^ = 1.2 ± 0.7 Rq and i?JF08 ^ 
12 ± 7 R0, respectively. The larger size of the progenitor 
of SN 2008D based on the model of Chevalier fc Frans- 
son is somewhat more consistent with typically observed 
radii of WN stars than the size required by the Waxman 
model. 

Furthermore, we find that the Swift X-ray spectrum 
can be fit equally well by an absorbed power law or a 
superposition of about equal parts of both power law 
and blackbody, which has not been suggested before in 
the literature, and has implications for determining the 
mechanism of the X-ray emission. Specifically, our de- 
rived blackbody radius fror n the PL-I-BB f i t is la rger (by 
a factor of 10) than that of iMazzali et al.l (|2008[ ). whose 
jet interpretation of the origin of the X-ray emission is 
based largely on the small inferred BB area. 

XRT 080109 is different from the apparent distribu- 
tion of classical X-ray fiashes as evidenced by its long 
duration, softening of the X-ray spectrum, and very high 
optical to X-ray luminosity ratio, though future studies 
of the inherent distribution will need to fold in detector 
effects and consider rest-frame properties. Future detec- 
tions are expected to show unambiguously whether these 
two events are part of a continuous distribution of XRFs, 
or members of a different physical class of events. 

Additional multi-epoch late-time observations (espe- 
cially for more than 200 days after maximum light) of 
SN 2008D will be useful to monitor the behavior of the 
double-peaked oxygen lines and their associated possible 
blueshift, and, as the SN becomes fully transparent, allow 
for closely constraining the geometry of the explosion as 
revealed by the core ejecta. Moreover, such observations 
are important to monitor for any late-time interactions 
between SN ejecta and the presumably prior-shed hydro- 
gen envelope to further study the mode of mass loss and 
search for episodes of sudden mass loss. 
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TABLE 1 
X-RAY Light Curve of XRT 080109 



t [s] 


dt [s] 


Luminosity Lx [10'*^ erg s ^] 


a^L^I [10^2 erg 


11.29 


11.28 


3.60 


1.06 


27.59 


5.02 


15.38 


4.53 


37.62 


5.02 


13.65 


4.34 


47.65 


5.02 


12.43 


4.14 


55.17 


2.51 


28.02 


8.67 


60.19 


2.51 


28.02 


8.67 


65.20 


2.51 


38.47 


10.11 


71.47 


3.76 


18.49 


5.78 


78.99 


3.76 


20.32 


6.04 


86.51 


3.76 


21.68 


6.29 


92.78 


2.51 


25.26 


8.28 


99.05 


3.76 


16.84 


5.52 


109.08 


6.27 


13.88 


3.92 


119.11 


3.76 


20.62 


6.04 


125.38 


2.51 


30.36 


9.06 


131.65 


3.76 


17.03 


5.50 


137.91 


2.51 


25.70 


8.23 


149.20 


8.78 


8.04 


2.48 


164.24 


6.27 


14.51 


3.91 


175.52 


5.01 


14.37 


4.34 


186.81 


6.27 


10.17 


3.30 


200.60 


7.52 


11.18 


3.14 


215.64 


7.52 


8.57 


2.75 


228.18 


5.01 


12.71 


4.13 


244.46 


11.28 


5.68 


1.83 


263.28 


7.52 


9.53 


2.89 


278.32 


7.52 


9.34 


2.89 


294.62 


8.78 


7.34 


2.36 


317.19 


13.79 


4.64 


1.57 


346.02 


15.04 


2.69 


0.86 


369.84 


8.78 


3.54 


1.18 


383.63 


5.01 


6.19 


2.06 


398.68 


10.03 


3.42 


1.09 


427.51 


18.81 


1.79 


0.58 


472.64 


26.33 


1.29 


0.41 


514.01 


15.04 


2.06 


0.69 


778.53 


249.47 


0.26 


0.060 


14944.11 


9184.50 


0.0040 


0.0055 


677205.64 


376223.75 


< 0.0016 




898641. 0** 


17900. 


0.00032 


0.00017 



Note. — Using a count-to-luminosity conversion of 4.3 X erg s""*^ per count s~^ 

based on tiie pure PL spectral fit, for a distance of D = 31 ± 2 Mpc. Swift data unless 
noted otherwise. 
^ Chanda data. 
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TABLE 2 

Model Parameters for X-ray Spectral Fits of XRT 080109 



Parameter Power-Law (PL) Fit Blackbody (BE) Fit Combined PL+BB Fit^ 



Quality of Fit {x^/dof) 


12.73/24 


28.56/25 


12.57/23 


Implied Nh [lO^^ cm-^] 


r 2+2.1 


< 1.5 


5 2+2'i 


Time-integrated F, kT, and [F;fcT] respectively 


9 1+0.3 
^•-^-0.4 


0.75 ± 0.07 


a.lloi; 0.10 ±0.01 


ix.iso*^ [ 10^2 erg s-i] 




4.6 ±0.6 


13 ± 2"= 


i^x.iso'' [10« erg] 


5.8l? 


2.4± 3 


6.7 ±4^= 



^ Assuming the same (factor unity) ratio of PL to BB flux as in XRF 060218 jButler et aU2006|) .^ Unabsorbed and timc-avcragcd 
value over the range 0.3—10.0 kcV (for PL fit) and bolometric (for BB fit), respectively, assuming — 31 ± 2 Mpc.*^ Listed only 
for the BB component. Lx,iBo and .Ex.iso are computed assuming — 10^^ cm as required from the combined PL and BB fit. 



TABLE 3 

Photometry of comparison stars in the field of SN 2008D by KAIT and FLWO 1.2 m 



ID 


U 


Nu 


B 


Nb 


V 


Nv 


R/r' 




I/i' 






[mag] 




[mag] 




[mag] 




[mag] 




[mag] 




1 


16.771(029) 


2 


16.438(002) 


3 


15.563(006) 


2 


15.117(013) 


3 


14.640(016) 


3 


2 


17.703(037) 


2 


17.732(003) 


2 


16.997(002) 


4 


16.590(005) 


2 


16.274(007) 


6 


3 


13.611(005) 


3 


13.745(008) 


3 


13.234(010) 


4 


12.933(010) 


4 


12.639(015) 


3 


4 


17.822(069) 


2 


17.865(009) 


3 


17.409(012) 


4 


17.061(010) 


5 


16.781(002) 


6 


5 






17.994(011) 


2 


17.316(010) 


4 


16.841(013) 


4 


16.498(009) 


5 


6 






18.132(014) 


4 


17.550(016) 


3 


17.115(005) 


4 


16.856(011) 


6 


7 


12.466(005) 


3 


12.530(002) 


3 














8 






18.977(012) 


4 


17.924(011) 


4 


17.227(006) 


5 


16.789(003) 


3 


9 


15.003(007) 


3 


14.969(003) 


3 


14.394(016) 


3 


14.041(010) 


3 


13.734(013) 


3 



1 




16.415(049) 


1 


15.586(033) 


1 


15.326(045) 


1 


15.088(046) 


1 


3 




13.657(031) 


1 


13.191(022) 


1 


13.089(031) 


1 


12.999(031) 


1 


10 




16.498(049) 


1 


15.594(033) 


1 


15.298(045) 


1 


15.047(045) 


1 


11 




16.908(051) 


1 


15.764(035) 


1 


15.344(047) 


1 


15.004(048) 


1 


12 




18.068(089) 


1 


17.237(060) 


1 


17.038(079) 


1 


16.805(081) 


1 


13 




15.659(040) 


1 


14.827(027) 


1 


14.585(037) 


1 


14.391(037) 


1 


14 




14.364(033) 


1 


13.631(023) 


1 


13.446(031) 


1 


13.311(031) 


1 


15 




17.584(086) 


1 


17.164(058) 


1 


17.087(076) 


1 


17.040(078) 


1 


16 




16.054(045) 


1 


15.320(031) 


1 


15.105(041) 


1 


14.946(042) 


1 


17 




17.309(055) 


1 


15.966(037) 


1 


15.486(049) 


1 


15.059(051) 


1 



Note. — Uncertainties (standard deviation of the mean) are indicated in parentheses. Stars (see Fig.[3j in the top 
half of the table were observed with KAIT in U BV RI and the stars in the bottom half were observed with the FLWO 
1.2-m telescope in BVr'i' . 
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TABLE 4 
Optical Photometry of SN 2008D 



MJD 






B 




V 


R 


fl(c)b 




r' 




i' 


/ 


Tcl.'^ 


[day 




[day] 


[mag 




[mag] 


[mag] 


[mag] 


[mag] 


[mag] 


[mag] 




54474 


45 


-0.12 


> 


20.10 


> 19.80 


— 


— 


> 


20.20 


> 


20.20 


— 


FLWO 


54475 


41 


0.84 


19 


08 


03) 


18.52(04) 


— 


— 


18 


31(04) 


18 


17(03) 


— 


FLWO 


54476 


41 


1.85 


19 


05 


07) 


18.33(04) 


17.80(03) 


17.78(05) 










17.34(05) 


KAIT 


54476 


43 


1.87 


18 


98 


03) 


18.38(03) 


— 


— 


18 


11(03) 


17 


86(03) 




FLWO 


54477 


35 


2.78 


19 


01 


03) 


18.39(05) 


17.84(03) 












17.36(02) 


Nickel 


54477 


36 


2.79 


19 


11 


02) 


18.44(03) 






18 


15(02) 


17 


86(02) 




FLWO 


54477 


46 


2.89 


19 


15 


08) 


18.43(05) 


17.85(04) 


17.82(05) 










17.38(04) 


KAIT 


54478 


34 


3.78 


19 


09 


03) 


18.37(02) 


17.89(02) 


— 










17.37(02) 


Nickel 


54478 


39 


3.83 


19 


24 


02) 


18.47(02) 


— 


— 


18 


17(02) 


17 


89(02) 


— 


FLWO 


54478 


44 


3.87 


19 


22 


09) 


18.43(05) 


17.83(04) 


17.81(05) 










17.40(04) 


KAIT 


54479 


43 


4.86 


19 


14 


06) 


18.30(04) 


17.78(04) 


17.76(05) 










17.19(04) 


KAIT 


54480 


39 


5.83 


19 


34 


06) 


18.22(04) 






18 


02(03) 


17 


72(04) 




FLWO 


54480 


40 


5.84 


19 


23 


07) 


18.22(04) 


17.68(03) 


17.64(05) 










17.17(03) 


KAIT 


54481 


43 


6.86 


18 


95 


11) 


17.97(05) 


17.47(06) 


17.48(05) 










16.97(05) 


KAIT 


54481 


51 


6.94 


19 


10 


04) 


18.05(04) 


— 




17 


75(05) 


17 


58(09) 


— 


FLWO 


54482 


39 


7.82 


18 


88 


1 1 ^ 
11) 


17.93(05) 


17.34(03) 


: 










16.81(04) 


KAIT 


54482 


43 


7.86 


18 


98 


02) 


17.94(03) 


— 


— 


17 


65(03) 


17 


34(03) 


— 


FLWO 


54483 


27 


8.71 


18 


89 


U5) 


17.70(03) 


— 


— 


17 


51(02) 


17 


25(03) 


— 


FLWO 


54483 


36 


8.79 


18 


90 


10) 


17.73(04) 


17.27(02) 


17.21(06) 










16.72(03) 


KAIT 


54484 


36 


9.79 


18 


69 


07) 


17.69(04) 


17.13(02) 


17.07(06) 










16.58(02) 


KAIT 


54485 


34 


10.77 


18 


54 


08) 


17.61(05) 


17.07(02) 


16.97(08) 










16.41(02) 


KAIT 


54489 


41 


19.90 


18 


76 


05) 


17.63(03) 


— 


— 


17 


00(02) 


16 


87(02) 


— 


FLWO 


54495 


36 


20.80 


18 


36 


04) 


17.38(03) 


— 


— 


17 


00(02) 


16 


65(02) 


— 


FLWO 


54496 


34 


21.77 


18 


52 


03) 


17.41(02) 


— 


— 


16 


95(02) 


16 


59(02) 


— 


FLWO 


54497 


39 


22.83 


18 


63 


03) 


17.40(02) 






17 


06(03) 


16 


68(03) 


— 


FLWO 


54498 


30 


23.74 


18 


63 


05) 


17.50(04) 






17 


05(02) 


16 


65(02) 




FLWO 


54499 


35 


24.79 


18 


72 


02) 


17.49(02) 






17 


07(02) 


16 


70(02) 




FLWO 


54501 


42 


26.85 


18 


92 


07) 


17.54(03) 


16.82(02) 












16.20(02) 


KAIT 


54502 


36 


27.79 


19 


01 


05) 


17.58(03) 


16.86(02) 












16.22(02) 


KAIT 


54502 


44 


27.87 


18 


93 


04) 


17.51(02) 






17 


18(02) 


16 


77(02) 




FLWO 


54503 


34 


28.78 


19 


15 


03) 


17.67(02) 






17 


22(02) 


16 


78(02) 




FLWO 


54503 


38 


28.81 


19 


13 


06) 


17.64(03) 


16.86(02) 












16.24(02) 


KAIT 


54504 


35 


29.78 


19 


19 


05) 


17.74(03) 


16.95(02) 












16.27(03) 


KAIT 


54504 


37 


29.81 


19 


25 


04) 


17.71(02) 






17 


29(02) 


16 


83(02) 




FLWO 


54505 


33 


30.76 


19 


28 


06) 


17.82(03) 


17.01(02) 












16.28(02) 


KAIT 


54507 


32 


32.75 


19 


52 


08) 


17.92(03) 


17.11(02) 












16.38(02) 


KAIT 


54509 


29 


34.72 


19 


70 


06) 


18.06(03) 


17.21(02) 












16.42(02) 


KAIT 



Note. — Uncertainties of the measurements arc indicated in parentheses. 
^ Time in days since to.^ i?(c)-band magnitudes derived from unfiltered observations.^ Telescope used: KAIT — 0.76- 
m Katzman Automatic Imaging Telescope; Nickel — Lick Observatory 1-m Nickel telescope; FLWO — FLWO 1.2-m 
telescope. 



TABLE 5 

Swift/VVOT PHOTOMETRY OF SN 2008D 



[day] 


<5tb 
[day] 


Exp " 
[s] 


Mag 


a 
[mag] 


Filter 


0.1359 


0.2688 


412.27 


20.24 


0.21 


U 


2.0047 


1.0058 


2677.07 


18.62 


0.03 


U 


2.9769 


0.8591 


1001.81 


19.07 


0.06 


U 


3.9479 


0.9349 


1314.54 


19.36 


0.07 


U 


4.9872 


0.9346 


1347.46 


19.37 


0.07 


U 


5.9550 


0.9326 


921.90 


19.45 


0.08 


U 


Note. — Table \E\ is published in its entirety in the 
eleetronic edition of the Astrophysical Journal. A por- 



tion is shown here for guidance regarding its form and 
content. Listed uncertainties do no include systematic 
uncertainties. 

^ Time in days (middle point of the combined image) 
since to^ 2008-01-09 13:32:49.^ Time bin in days {from 
the start of the first exposure to the end of the last ex- 
posure in the combined image). Total exposure time (in 
seconds). 
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TABLE 6 

Near-Infrared Photometry of SN 2008D 



MJD 




Ac 




J 


— 

H 




Kg 


r J 1 
[day] 




r J 1 
[day] 


[mag] 


[mag] 


[mag] 


54474, 


35 


-0.21 


> 


18.15 


> 17.12 


> 16.18 


54475. 


27 


0.71 


17, 


.00(03) 


16.70(05) 


17.09(12) 


54476. 


33 


1.77 


16, 


.54(07) 


16.22(09) 


16.22(10) 


54477. 


,34 


2.78 


16, 


.45(07) 


16.15(09) 


16.13(15) 


54478. 


,34 


3.78 


16, 


.44(04) 


16.05(09) 


15.99(12) 


54479. 


,38 


4.82 


16, 


.38(10) 


15.96(10) 


15.81(10) 


54480. 


,34 


5.78 


16, 


.23(10) 


15.86(10) 


15.63(07) 


54481. 


,32 


6.76 


15, 


.99(10) 


15.72(07) 


15.55(02) 


54482. 


,26 


7.70 


15, 


.88(04) 


15.55(00) 


15.46(05) 


54483. 


,32 


8.76 


15, 


.74(08) 


15.41(13) 


15.21(06) 


54484. 


,35 


9.79 


15, 


.62(05) 


15.31(05) 


15.11(07) 


54485. 


,31 


10.75 


15, 


.51(02) 


15.17(16) 


15.05(13) 


54486. 


,35 


11.79 


15 


.43(01) 


15.12(05) 


14.92(05) 


54497. 


,19 


22.63 


15 


.00(04) 


14.76(08) 


14.53(03) 


54499. 


,26 


24.70 


15 


.07(21) 


14.73(03) 


14.44(04) 


54505. 


,20 


30.64 


15. 


.11(05) 


14.76(02) 


14.54(11) 


54508. 


,24 


33.68 


15. 


.25(02) 


14.86(01) 


14.68(13) 


54509. 


,22 


34.66 


15. 


.24(05) 


14.88(03) 


14.68(05) 


54510. 


,25 


35.69 


15. 


.32(04) 


14.89(06) 


14.68(03) 



TABLE 7 

Journal of Spectroscopic Observations of SN 2008D 



UT Date 


At"" 


. b 
••Vmax 


Telescope*^ 


Range'* 


Airmass° 


Slit 


Exp. 








[days] 


[days] 




[A] 




["] 


[sec] 


2008 


Jan. 


11.27 


1.70 


-16.7 


GMOS-S 


4000-8500 


2.25 


0.75 


3 x600 


2008 


Jan. 


11.41 


1.84 


-16.6 


MMT 


3200-8400 


1.12 


1.00 


1200 


2008 


Jan. 


12.36 


2.80 


-15.6 


Lick 


3300-8100 


1.04 


2.00 


2x 1200 


2008 


Jan. 


12.55 


2.99 


-15.4 


MMT 


3300-8100 


1.62 


1.00 


2x 300 


2008 


Jan. 


12.67 


3.11 


-15.3 


Keck II 


5000-10000 


1.64 


1.00 


300 


2008 


Jan. 


13.30 


3.74 


-14.7 


MMT 


3300-8100 


1.08 


1.00 


1200 


2008 


Jan. 


13.38 


3.82 


-14.6 


Lick 


3300-8100 


1.01 


2.00 


2x 1200 


2008 


Jan. 


14.36 


4.79 


-13.6 


MMT 


3300-8400 


1.00 


2.00 


1200 


2008 


Jan. 


14.40 


4.84 


-13.6 


Lick 


3300-8100 


1.01 


2.00 


4x 1200 


2008 


Jan. 


14.45 


4.89 


-13.6 


IRTF 


8000-25,000 


1.07 


0.5 


10 x 150 


2008 


Jan. 


15.31 


5.75 


-12.7 


MMT 


3300-8400 


1.05 


3.50 


1800 


2008 


Jan. 


15.37 


5.80 


-12.6 


Lick 


3300-10700 


1.02 


2.00 


2x 1800 


2008 


Jan. 


16.34 


6.77 


-11.7 


MMT 


3300-6500 


1.05 


1.00 


3x 900 


2008 


Jan. 


18.57 


9.00 


-9.4 


Lick 


3500-10700 


1.80 


2.00 


2x 1800 


2008 


Jan. 


20.27 


10.70 


-7.7 


APO 


3500-5600 


1.13 


1.50 


2x 1200 


2008 


Jan. 


21.22 


11.66 


-6.8 


APO 


3300-9800 


1.30 


1.50 


3x 1200 


2008 


Jan. 


21.51 


11.95 


-6.5 


IRTF 


8000-25,000 


1.04 


0.5 


11 X 150 


2008 


Jan. 


28.51 


18.95 


0.5 


IRTF 


8000-25,000 


1.07 


0.5 


11 x 150 


2008 


Feb. 


01.35 


22.78 


4.3 


FLWO 


3500- 7400 


1.00 


3.00 


2x 1800 


2008 


Feb. 


02.44 


23.88 


5.4 


APO 


3300-9800 


1.25 


1.50 


2x 1200 


2008 


Feb. 


09.12 


30.56 


12.1 


MMT 


3300-8400 


1.05 


1.00 


3x 900 


2008 


Feb. 


11.18 


32.61 


14.2 


MMT 


3300-8400 


1.05 


1.00 


3x 900 


2008 


Feb. 


12.43 


33.89 


15.4 


Keck 


3100-9300 


1.03 


1.00 


2x 900 


2008 


Feb. 


16.36 


37.79 


19.4 


Lick 


3300-10700 


1.02 


2.00 


1800 


2008 


Feb. 


27.08 


48.52 


30.1 


LOO 


3800-10800 


2.57 


1.0 


3x300 


2008 


Feb. 


28.37 


49.81 


31.4 


APO 


3300-9800 


1.24 


1.50 


2x 1200 


2008 


Mar. 


10.23 


60.67 


42.2 


Keck 


3100-9100 


1.30 


1.00 


600 


2008 


Mar. 


14.48 


63.92 


45.5 


IRTF 


8000-25,000 


1.54 


0.5 


5 x 150 


2008 


Apr. 28.33 


109.77 


91.4 


Keck 


3100-9300 


1.30 


1.00 


2x 400 


2008 


Jun. 


07.26 


149.69 


132.4 


Keck 


3100-9300 


1.50 


1.00 


2x 900 


Note. — 


This Tabic is published in its entirety in tlie electronic edition of the Astrophysical Journal. 



A portion is shown here for guidance regarding its form and content. 

^ Days after X-ray outburst, to^ 2008-01-09 13:32:49 UT.^ Days with respect to date of V- 
band maximum.'^ APO — 3.5m Telescope at Apache Point Observatory/DIS; FLWO — Tillinghast 
1.5m/FAST; GMOS-S = Gemini-South 8.4m/GMOS; IRTF^ IRTF 3m/SpeX; Keck I = Keck I 
lOm/LRIS; Keck II = Keck II lOm/DEIMOS; LGO = 6.5m Clay Magellan Telescope at Las Gam- 
panas Observatory/LDSS-3; Lick = Lick 3m/Tillinghast; MMT ^ 6.5m MMT/Bluechannel.^ Airmass 
at the middle of the set of observations.® Observed wavelength range of spectrum. In some cases the 
blue and red ends of the spectrum are extremely noisy and are not shown in the figures. 
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TABLE 8 

Blackbody Fits and Bolometric Light Curves of SN 2008D 
PROM UVWl -> Ks Data 









Rbb 


log{LBB) 


log(L| 


UVWl-K.) 


[davsl 


fKl 


[IQl-l cm] 


[ erg s-i] 


[ erg s 1] 


0.80 


12000l-«» 


9 yt-+0.36 


42, 


,04+°-2J 


41, 


09+0. 18 
17 


1.80 


lOOOOtf™ 


o vo+0.30 
■^■'^-0.30 


42, 


,05+f,:« 


42, 


-oo+°-l^ 


2.80 


8700llg«« 




41, 


,93+;;:!^ 


41, 




3.80 


7900110°" 


5. 


-10+0.47 
^- -0.47 


41, 




41, 


00 + 0. 14 
13 


4.80 


7600l?^o 




41, 




41, 


85+°-14 
12 


5.80 


7400+^1° 


6. 


19+0.56 
^^-0.54 


41, 


gQ + 0:i4 

— 12 


41, 


og+o:i3 

~0 12 


6.80 




g 00+0.61 


41, 


g5+o:i3 


41, 


93+0:13 


7.80 


7000l»'» 


7 64+°-65 
'•"^-0.64 


42, 


-oo+°-i! 


41, 




8.80 


69001^™ 


»-3U_o gg 


42, 


,04+°:i? 


42, 


.02+°:i2 


9.80 


6800l^to 


8.89_q;73 


42, 


,08+°:i2 


42, 


•05+on 


10.80 


67001^20 


''•^'-'-0.76 


42, 




42, 


.08+°-i? 


11.80 


66001™ 


•-'-(). 81 


42, 


1 n + 0.12 
Oil 


42, 


1-1+0:12 

11 


12.80 


6600+^f(! 


10 


40+0-84 


42, 


1 c+o:i2 

() 10 


42 


-1 0+0:12 
{) 11 


13.80 


6500l««° 


10 


sn+0-92 

•*U_o.84 


42, 




42 


-i i-+o:i2 

Oil 


14.80 


6400l«i0 


11, 


90+0-94 
■"^"-0.86 


42, 


-is+o-ii 

■lO Q -J^Q 


42 


17+0:12 

Oil 


15.80 


6400l«30 


11, 


gQ+0.90 


42, 


19 + 0:11 

10 


42 


•18+on 


16.80 


6300l«2O 


11 


90+0.97 


42, 


20+°-ll 
10 


42 


19+0:11 

Oil 


17.80 


6200+^5'(J 


12 


on+0.98 
•■^^-0.96 


42, 


91+0:11 

J- 


42 


20+0:11 

10 


18.80 


620011^° 


12 


60+1-"" 


42, 


91+0:11 

-^-'- — 0.10 


42, 


20+"-ii 

-^'J — 0.10 


19.80 


eiooli™ 


12, 


00+1-°° 


42, 


91+0.11 
■^^-0.10 


42, 


20+°-ii 
•^'-'-0.10 


20.80 


eooollfg 


13 


9n+i-00 

•^'^-0.99 


42, 


91+0.10 
■^J^-0.10 


42 


2Q+0.11 
•^'-'-0.10 


21.80 


6000l2« 


13 


40+1-1" 


42, 


9O+0.10 
■^^-0.09 


42 


20+0.11 

•^'^-o.io 


22.80 


5900lif(; 


13 


7n+i.()() 

•"^-1.00 


42, 


iq+0.10 

■-^-'-0.09 


42 


-19+0-11 

--"-■-•-o.io 


23.80 


58001^^0 


13 


qn+i.io 


42, 


10+0.10 

■-^"-0.09 


42, 


iQ + O-lO 

--■-o-o.io 


24.80 


570011°° 


14, 


20+11° 


42, 


17+0.10 
■^'-0.09 


42, 


16+°-l° 

•J^°-o.io 


25.80 


5600l«° 


14, 


50+1-1° 


42, 


ir,+0.09 
■^"^-0.09 


42 


-TC+O-IO 
•^"^-0.09 


26.80 


5400l«[! 


14 


00+1.20 


42, 


1 O+0.09 
■-^"'-0.08 


42 


1 o+O.lO 

•^■'-o.og 


27.80 


5300l^«^ 


15 


90+1.20 
•^^-1.20 


42, 


11+0.09 
■^^-0.08 


42 


1 1+0.10 
•^^-0.09 


28.80 


5200l«° 


15 


50+1-2" 


42, 


Qg + 0.09 

■"'^-0.08 


42 


09+0-10 

•"-^-0.09 


29.80 


5100l«° 


15, 


■'"-1.20 


42, 


00+0. 09 
■U»_0.08 


42, 


Oo+O-lO 
•'J**-0.09 


30.80 


5100l«° 


15, 


O0+1-20 
•°U_i.20 


42, 


Qg+0.08 
■'J°-0.07 


42 


07+0-09 
•"'-0.09 



Note. — Adopted values: E(B - V')HoBt = 0.6 ± 0.1 mag, D = 31 ± 
2 Mpc, and to= 2008-01-09 13:32:49 UT. The uncertainties include the 
uncertainty in reddening, but not the systematic uncertainty in distance. 



TABLE 9 

Optical and Near-Infrared He I Absorption-Line Velocities of SN 2008D 



JD-2,400,000 


At 


^Vmax 


He I A5876 


He I A6678 


He I A7065 


He I 1.0830 /im 


He I 2.0581 /um 


[day] 


[day] 


[day] 


[km s-l] 


[km s-l] 


[km s-l] 


[km s-l] 


[km s-l] 


54479.96 


4.89 


-13.6 








-21,700^ 


-14,200 


54480.82 


5.75 


-12.7 


-14200 


-11800 


-10500 






54481.84 


6.77 


-11.7 


-13300 


-12700 


-13000 






54484.07 


9.00 


-9.4 


-12300 


-12400 


-10900 






54486.72 


11.66 


-6.8 


-11500 


-11100 


-10500 






54487.01 


11.95 


-6.5 








-14,900 


-12,200 


54494.01 


18.95 


0.5 








-13,900 


-11,700 


54497.85 


22.78 


4.3 


-10600 


-10400 


-10100 






54498.94 


23.88 


5.4 


-10600 


-10200 


-10200 






54505.62 


30.56 


12.1 


-10500 


-10200 


-9900 






54507.68 


32.61 


14.2 


-10400 


-10000 


-9900 






54508.93 


33.87 


15.4 


-10300 


-9900 


-9700 






54512.86 


37.79 


19.4 


-10000 


-9600 


-9500 






54524.87 


49.81 


31.4 


-9300 


-8600 


-8700 






54535.73 


60.67 


42.2 


-8600 


-7900 


-8100 






54539.0 


63.92 


45.5 








-13900 


-10200 



^ Likely due to a blend and not traced by He I. 



